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RATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2130 

CALCULATION OF TRANSONIC FLOWS PAST THIN 
AIRFOILS BY AN INTEGRAL METHOD 
By William Perl 


■SUMMARY 

A method of calculating two-dimensional compressible flows past 
thin airfoils is presented with particular reference to the transonic 
speed range. The method is based on the integral form of the equa- 
tion of continuity and on the intrinsic form, in terms of the stream- 
line curvature, of the irrotationality condition. The application to 
the transonic range conforms with, and is carried out within the 
framework of, the transonic similarity theory. The results include 
velocity distributions for various symmetric sections in the 
continuous-potential-flow speed range and for the symmetric biconvex 
airfoil at zero angle of attack through the transonic speed range 
and the associated variation of terminal- and head-shock location 
and pressure drag with Mach number. Comparisons with other theories 
and with available experimental data indicate at least qualitatively 
good agreement. The location and the movement of the terminal shock 
in the lower transonic speed range are discussed, as well as the 
limiting- line phenomenon. 


INTRODUCTION 

Considerable effort has been expended in recent years on the 
investigation of the flow phenomena that occur when aerodynamic 
bodies travel at transonic speeds. Experimentally, the main fea- 
tures of steady transonic flow patterns about airfoil sections have 
been established. (See, for example, references 1 to 5.) 

Thus, the succession of flow patterns, of which an approx- 
imate theoretical description is attempted herein, may be briefly 
outlined as follows : Consider the zero -lift flow relative to . 

an airfoil that is symmetric with respect to its chord line 
(fig. l) . At low subsonic free-stream velocities, the flow 
pattern is likewise symmetric with respect to the chord line and is 
everywhere subsonic. Above a certain subsonic free-stream velocity, 
a symmetric-type local supersonic region, bounded by the locus of 
sonic speed indicated by 1 (fig. 1), starts to develop in the 
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neighborhood of the maximum- thickness location of the airfoil. If 
the airfoil had fore-and-aft symmetry with respect to its midchord 
location, these local supersonic regions would likewise have fore- 
and-aft symmetry with respect to midchord. As the free- stream 
velocity is further increased, the fore-and-aft type of symmetry 
of the flow pattern disappears and an asymmetric pattern develops. 

The local supersonic region, indicated by the sonic boundary 2 in 
figure 1, is now of an asymmetric type (that is, would be asymmetric 
fore and aft even if the airfoil were symmetric fore and aft) and 
is terminated on the downstream side by ah approximately normal 
shock. Hear the airfoil especially, this terminal shock is actually 
a more or less complex shock structure depending on the Reynolds 
number and interacts with the boundary layer in a complex manner. 

As the free-stream velocity increases toward sonic, the asymmetric- 
type local supersonic region increases in extent and the terminal 
shock' increases in intensity and moves downstream toward the trail- 
ing edge. At slightly supersonic free-stream velocities, an addi- 
tional shock, the bow wave, appears ahead of the airfoil and 
approaches it as the free-stream speed is increased. The subsonic 
part of the flow field has now become finite in extent, as indicated 
by the sonic boundary 3 (fig. 1) . The terminal shock has reached, 
and remains in, the neighborhood of the trailing edge and begins to 
slope in a downstream direction. At a certain supersonic free-stream 
6peed the forward shock attaches at the nose of the airfoil, if the 
nose is sharp-edged, and shortly thereafter the flow field becomes, 
usually, everywhere supersonic, except, of course, in the boundary 
layer. 


On the theoretical side, much of the effort expended on the 
transonic problem has gone into the solution of the mathematically 
well-defined problem of the continuous compressible potential 
flow past an airfoil, corresponding to the symmetric- type flow 
pattern described. (See, for example, references 6 to 20.) Emmons 
and Maccoll and Codd (references 18 to 20), however, have included 
shocks in some of their relaxation calculations. Progress with the 
methods cited has so far not been rapid, primarily because of mathe- 
matical and computational difficulties. Perhaps the most significant 
theoretical result that lias emerged is the probability that a sub- 
sonic free-stream speed exists for a given aerodynamic shape, such 
that continuous symmetric-type flow patterns containing local super- 
sonic regions can, at least in principle, be derived for this shape 
below this speed but not above. In regard to the other main 
ingredient of real transonic flows, namely, the boundary layer and 
its interaction with the terminal shock. Interpretation and analysis 
is only in a prelim! nary stage. (See, for example, references 3, 

4, and 21 to 24.) It must be noted, however, that progress with a 
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complete theory is hampered hy the practical difficulties involved 
in carrying out good experiments in the transonic speed range. 

In the present report an integral method is applied to the cal- 
culation of transonic flow patterns. Starting from the equation of 
continuity in integral form and conservation of energy in the inte- 
grated form of Bernoulli’s equation, which suffice for the analysis 
of quasi one-dimensional flows, two-dimensional effects are then 
taken into account hy inclusion of the equation of motion normal to 
the streamlines, expressed in terms of the streamline curvature . 

By an appropriate choice of the streamline curvature, the flow 
pattern can then he calculated in a comparatively simple manner. 

This method has long been known and applied to the calculation 
of flows through curved channels, for example, between turbine 
blades, or in wind tunnels (references 25 and 26). Application has 
recently been made to isolated airfoils in continuous compressible 
potential flow (references 27 to 29). 

The small -perturbation form of the method of reference 29 is 
developed herein. The resulting simplifications make feasible the 
approximate calculation of some of the more complicated, shock- 
containing, flow patterns previously described. The small -perturbation 
form also shows explicitly the agreement of the method of reference 29 
with the more general transonic similarity theory of von Barman 
(reference 30) and Oswatitsch, to which any ncnviscous-type explicit 
solution for thin airfoils in the transonic speed range must conform. 

The mathematical simplicity of the present method preserves the 
possibility of interpreting results and difficulties at all stages 
in terms of the basic conservation laws. The insight thus permitted 
into the nature of transonic flows makes it easy to understand, for 
example, the previously mentioned phenomenon of an upper- limiting 
subsonic flight speed for the existence of continuous potential 
flows . 

In order not to unduly disturb the continuity of presenta- 
tion of the main ideas, certain mathematical details Eire relegated 
to appendices A to D and a list of. the main symbols used is con- 
tained in appendix E. 

The main developments of this report were worked out at 
Columbia University in the winter of 1947-1948 under the stimulus 
of a lecture given by Professor Theordore von Karman on the subject 
of transonic similarity. The report was completed at the NACA Lewis 
laboratory o 
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The author wishes to express his thanks to Professor Theodore 
von KvCrman for having inspired the present treatment of the tran- 
sonic problem. This inspiration was not only direct; for, Professor 
von Karman has probably done most to demonstrate the power and 
utility in fluid mechanics of the two foundational ideas on which 
the present theory is built , namely, the integral-method approach 
and the similarity concept. It is to be hoped that the present 
structure does not rest too lightly on the foundations. 


COMPRESSIBLE POTENTIAL FLOW 

Governing equations. - The flow is considered steady, two- 
dimensional, isentropic, irrotational, and symmetric with respect 
to the airfoil chord (x-axis). The equations of motion can be 
expressed in terms of the local fluid density p, speed v, pres- 
sure p, Mach number M, the corresponding free- stream values 
(denoted by the subscript o), the specific-heat ratio 7 , and 
the streamline curvature C, as follows (fig. 2): (The main symbols 

are defined in appendix E.) 

Equation of continuity: 


rjc 


Pn 


PoVar = 


pvdn 


0 


Irrotationality condition (reference 31) : 

+ Cv = ° 

Bernoulli- state equation: 


( 1 ) 


( 2 ) 



( 3 ) 


The differentials dy and dn in equation (l) are elements of 
length along a potential line in the free stream far ahead of the 
airfoil and along a potential line in the region of the airfoil, 
respectively. The lower limits of integration lie on the zero 
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streamline, that is, on the s-azis and on the airfoil contour. The 
upper. limits lie on a streamline .ECEEF, ■which, as y Q approaches 

infinity, becomes parallel to the x- axis'. It is assumed for the 
present that as y 0 approaches infinity the displacement component 
5 in the y-direction between any two points on the same streamline, 
such as D and B, approaches zero. This assumption is valid for 
circulation- free (zero-lift) flows that are subsonic at infinity 
and that have no wakes or unclosed portions of airfoil contour 
extending to infinity in the x-direction, inasmuch as such flows 
can be represented at infinity by a collection of incompressible- 
type source- sink doublets in the finite part of the plane, and for 
such doublets the assumption is valid. This assumption can also be 
regarded, perhaps more generally, as restricting the flow pattern, 
to that produced in a closed, straight-walled, two-dimensional wind 
tunnel in the limiting case in which the walls are infinitely for 
apart. 


The element of length dn in equation (l) is now approximated 
by the element of length dy in the y-direction, that is, the 
potential lines are replaced by the lines x = constant. This 
replacement is exactly valid for the potential line starting from 
the midchord location o*f a symmetric airfoil that has also fore- 
and-aft symmetry. In appendix A, this replacement is shown to be 
justified, in general, . for the case of small perturbation of the f ree- 
stream flow by an airfoil. 


Next let the lower limit on the right side of equation (l) start 
from a point on the airfoil at which the airfoil ordinate is T. 

Set p Q , v Q , and p Q equal to unity, so that the corresponding 

local quantities p, v, and p henceforth represent nondimensional 
fractions of free- stream density, velocity, and pressure, respec- 
tively. Furthermore, let the unit of length be the airfoil chord c 
6o that the lengths x, y, and l/C henceforth represent non- 
dimensional fractions of the airfoil chord length. Subtracting 
y 0 -Y from both sides of equation (1), letting y Q approach infin- 
ity, and using the preceding assumption for 6 and approximation 
for dn yields 

poo 

T = / (pv-l)dy (4) 

Equations (2) and (3) become, respectively. 


dy = - 


dv 

Cv 


(5) 
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The form of the continuity equation (4) is analogous to that defining 
the displacement thickness in the "boundary-layer theory. The airfoil 
ordinate T here takes the place of the "boundary- layer displacement 
thickness. Whereas the subtraction and limiting process leading to 
equation (4) is performed in the boundary- layer theory mainly for 
definiteness in defining a boundary- layer thickness, the process is 
here essential in order to subtract out the infinite mass flow 
crossing a potential line and initially introduced perhaps somewhat 
artificially by considering the fluid as moving relative to the 
stationary airfoil. 

The streamline curvature C, considered positive when the 
streamline curves in the same sense as the corresponding contour of 
the airfoil at the maximum thickness location, is a function, of the 
position coordinates x and y of the flow field. At a given 
chordwise location x, the streamline curvature is a function of y 
that varies continuously from the curvature C a of the airfoil 
contour at y = T to the value zero at y = a>. It proves conven- 
ient to consider the curvature C at a given x as a function of 
the local velocity v rather than of y. In the applications made 
herein, v will be a monotonic function of y. Hence, at least for 
such applications, C can be expressed as 

C = C a f(x,v) " (7) 

The assumption that the streamline curvature can be as well 
expressed by a function of x and v as by a function of x and 
y may be more generally justified as follows: Of the four basic 

variables of the problem, namely, the vector velocity given by the 
two variables v and 9 , where 9 is the angle of Inclination of 
the local stream direction to the x-axis, and the vector distance 
given by the two variables x and y, two variables are to be 
regarded as independent and two as dependent. The solution of the 
so-called direct problem of airfoil theory is required in the form 
where v(x,y), 0(x,y), and in particular the velocity Y(x,T) on 

the given airfoil contour 0 a (x,Y) are finally known. ’ The analysis 
of the problem, however, can be made, at least in principle, with 
any two of the variables considered as independent and the remaining 
two as dependent. The hodograph method, for example, is based on v 
and 9 as the independent and x and y as the dependent variables. 
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The present method is based on v and x as the independent and 8 
and y as the dependent variables. Conceivably, other possible 
combinations might be advantageous in particular problems. Similar 
remarks apply with the potential, or the stream function, or both, 
replacing one or two of the present variables. 

The present method is thus a sort of semihodograph method. 
Whereas it retains the principal advantage of the hodograph method, 
namely, that the velocity v and hence, by equation (6), the density 
p becomes essentially an independent variable, it avoids the well- 
known disadvantage of the hodograph method; namely, difficulty in 
satisfying the boundary conditions in obtaining flows past given 
fixed airfoils as a function of Mach number. 


( 7 ) 


y 


Equation (4) becomes, upon substitution of equations 


(5) and 


YC a 




1 

p " v 
f(x,v) 


dv 


( 8 ) 


in which Y is the sought-for velocity at the airfoil at chordwise 
location x, the quantity x being held constant in the integration. 
Equations (6) and (8) constitute an expression for the velocity Y 
at a point on an airfoil in terms of the shape parameter YC a of the 
airfoil at that point, the free- stream Mach number Mq, and the 

curvature function f(x,v). (From the point of view of the inverse 
problem, that is, designing an airfoil to satisfy a prescribed 
velocity distribution, the shape parameter YC Q would be regarded 

as dependent upon Y, Mq, and f.) The basic reason for expressing 

the airfoil velocity Y in terms of a streamline- curvature function 
f is that, in the form of equation (8), Y is relatively insensitive 
to choice of f . As will be seen, a reasonable choice of f can 

therefore be made that leads with comparatively little computation to 

fairly accurate values of Y. 


The relation between the local velocity v and the lateral 
distance y at a given x is obtained by integration of equation (5) 
from an airfoil point (x,Y) or (x,Y) to a field point (x,y) or 
(x,v) as 


y = * + tt 

u a 


OY 


dv 

rf(x,v) 


( 9 ) 
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Solution of equation (8) for V enables the local velocity v(x,y) . 
at any point in the flow field to be determined from equation (9). 

The streamline inclination 0(x,y) can then also be obtained, as 
will be subsequently shown. 

The only approximation thus far made is to replace the element 
of length dn along a potential line by the element of length dy 
perpendicular to the chord line. The resulting equations (6), (8), 
and (9), together with an assumed form for the curvature function 
f(x,v), have been applied in reference 29 to calculate (finite 
perturbation) compressible potential flows past airfoils. An 
approximation for the density p(v) is next made, which leads to 
the general small-perturbation form of the preceding theory. 

Small -perturbation case. - By small perturbation is herein 
meant that the local velocity increment v-1 produced by the airfoil 
is small relative to unity. The density p is, by equation (6), 
thereby expressible as a power series in v-1. Thus, the numerator 
of the integrand in equation (8) can be developed as 

p - i * HV - r_u 2 + . . . (10) 

v c 


where 


V av - 1 


(11a) 


^l-«o 2 




(lib) 

(lie) 


iim r c s r = 2— did) 

Mo-»l 


The velocity increment A at the airfoil, a normalized local- 
velocity variable z, and a Mach number parameter K are next 
defined as 
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A = V - 1 


_ V T - 1 

A "V - 1 


l-M, 


K = 


(jT) 2 ^ (tP) 


o 

2J5 


(12a) 

(12b) 

(12c) 


In terms of the reduced local velocity z, the curvature function 
(equation (7)) can be expressed as 

C = C a f(x,z,K) (13) 

Equation (13) can easily be shown, from the general transonic sim- 
ilarity theory of reference 30, to be the most general possible 
form for the streamline- curvature function in the small-perturbation 
transonic range. (This possibility was pointed out by Professor 
Theodore von Kerman. ) Conversely, it will be presently verified 
that the use of equation (13) with the present formulation of the 
equations of motion leads to the transonic similarity laws of refer- 
ence 30. The curvature function f(x z,K) in equation (13) 
must satisfy the boundary conditions f(x,l,K) = 1 at the airfoil 
and f(x,0,K) =0 at infinity. 


Substitution of equation (10) Into equation (8) and changing 
the variable of integration from v to z by equations (lla), 
(12), and (13) yields 


PI 



£cA z 2\ dz 
p ) f(x,z,K) 


(14) 


Similarly, the relation (9) between local velocity v and lateral 
distance y becomes, for small perturbations. 


y 



J z 


dz 

f(ac,z,K) 


(15) 


Prandtl-Glauert range, - In the subsonic or Prandtl-Glauert 
range, defined by 0 < Mq < < 1 or E-> 00 , the dependence of the 
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curvature function (equation (13)) on K must disappear. Hence, 
in the small-perturbation case, A < < 1 , and for sufficiently 
subsonic free- stream Mach numbers M Q , the second term in the 
integrand of equation (14) can be- neglected relative to the first, 
and equation (14) can be written as 

YC„ 

— | = ?(x) (16) 

mA “ 


where F(x) is a function of x of the order of magnitude of 
unity. For smoothly curved, thin airfoils, the shape parameter YC a 
can be expressed as the product of the square of the thickness ratio 
T (T = maximum thickness /chord length) and a function of x of the 
order of magnitude of unity over most of the airfoil. Equation (16) 
can therefore be written as 

A = -~z_ G(x) (17) 

Vd 

where T << 1 and G(x) ~ 1. Equation (17) expresses the basic 
result of thin-airfoil theory that the velocity increment at chord- 
wise location x of one of a series of thin airfoils, all having 
the same geometric shape but differing in thickness ratio t, is 
proportional to the thickness ratio r of that airfoil. Equa- 
tion (17) also expresses the Frandtl-Glauert rule that the effect 
of compressibility in the lower subsonic rang e is to increase the 

airfoil velocity increment by the factor 1/V 1-Mo 2 , 

Lower transonic range; similarity rules. - As the free- stream 
Mach number Mq is increased toward unity, or p — > 0 , the airfoil 
contour being such as to maintain a small velocity perturbation 
A < < 1, a value of p < < 1 must eventually be reached where 
TA/p becomes comparable to unity. The second term in the integrand 
of equation (14) is then comparable with the first. Additional 
terms in the expansion of the integrand, however, evidently would 
remain negligible relative to the first two terms. Hence, at tran- 
sonic speeds defined by p ~ 0, instead of the relation of equa- 
tion (16) between the two parameters x and YC a /pA 2 , equa- 
tion (14) constitutes a relation among four parameters; thus 


H 


(18) 
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The last three parameters indicated in equation (18) are not inde- 
pendent, however, hut can he expressed in terms of two independent 
parameters. Thus, using the third parameter of equation (18) to 
eliminate A from the second yields the equivalent expression of 
equation (14) 


I 



u rA\ 

( T r) 2 7 3 ' p.J 


= o 


(19) 


or also, using the second parameter of equation (19) to alter the 
third. 


J 


I273 . ) 


= 0 


( 20 ) 


Inasmuch as the pressure coefficient P at a point on the airfoil 
is obtained in terms of the velocity from equation (3) in the small- 
perturbation case as 


Pi , P ~ P - ° p « -2 A (21) 

I Po-o 

relation (20) can he written as 

t 2 / 3 , * 

p " ^ M (22) 

The function in equation (22) is determined by the shape of the 
airfoil and is of the order of magnitude of unity. Equations (22) 
and (12c) agree with the transonic-similarity theory as regards 
velocity distributions on airfoils* 

By multiplying equation (15) by and using the preceding 

relations, there results for the velocity increment in the flow 
field 


= Jl^i&' ^ ( y- T) > (Tl) l/3 (y-T), k] (23) 


fo73 


V 
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where and *^2 are functions of the order of magnitude of unity 

and a/^T( 3 t— x) , the affinely transformed lateral distance of the 
Barman theory (reference 30) is of the order of magnitude of unity 
by equation (15). Equation (23) agrees with the transonic-similarity 
theory as regards conditions in the flow field. Equation (23) 
reduces to equation (22) on the airfoil. 

Equation (22) states that if the pressure coefficient has a 
certain value P^_ at a chordwise location z on an airfoil of 

thickness ratio in a flow at speed M Q ~ 1 and in a fluid 

of specific-heat ratio 7^, then the pressure coefficient Pg at 

the same chordwise location of a geometrically similar airfoil of 
thickness ratio Tg in a fluid of specific-heat ratio 7 2 and at 

a speed M Q ^ 2 given by 


will be related to P^ by 




Pi 



2/3 / r a/3 

' yv 


(25) 


The Prandtl-Glauert rule is included in equation (22) if an inverse 
square root is chosen for the K-dependent part of the function 3 * . 
This result necessarily follows because the Prandtl-Glauert case is 
obtainable from the general formulation, equations (14) and (19), 
by letting TA/d 0. 

A further similarity rule that concerns the local Mach number 
may be noted. In isentropic flow, the local Mach number M is 
related to the local velocity v and the free- stream Mach number 
Mq as follows: 

M V 

m 2 . 2 (26) 

1 - M 0 2 (v 2 -l) 

In the small-perturbation case v-l= v << 1, equation (26) 
becomes 
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M 2 -1 = - \i + 2T m V (27) 

where 

r M “ “o 2 (* + HT M o 2 ) ■ < 28a > 

and 

, lim r M = r = 2±1 ( 28 b) 

Mo-*! 

Equations (23), (27), and (28b) yield the transonic similarity rule 
for local Mach number 


(Tr) X / 3 (y-Y), k] 


(rr) 


(29) 


Equation (29) may be applied in a manner similar to equations (24) 
and (25) to relate local Mach numbers M-j_ and Mg at the 

corresponding points 


= Xg (30a) 

(T 1 r i) l/3 (y 1 -T 1 ) = (T^) 1 /*^) (30b) 

in the flow fields of similar airfoils and Tg in fluids 

characterized by and Tg at corresponding free-stream Mach 

numbers Mq ^ and Mg g given by equation (24). Instead of equa- 
tion (25), the following relation, obtained from equation (29), 
applies: 


Mg 2 -1 % 2 -l 

(T 2 r 2 )2/3 " ( T ! r i) 2/3 


(31) 


Various special cases of this similarity rule have been previously 
given in references 4, 32, and 33. 
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It may be noted that although the use of T= (7+l)/2 in such 
expressions as equations (23) and (29) is valid in the transonic 
limit Mq = 1, in applications for which Mq is finitely different 

from unity equations (28a) and (lie) for and r c , respectively, 

should more properly be used. The reason is that in the approxi- 
mating equations (10) and (27), only the small -perturbation approxi- 
mation V<< 1 was made and not the transonic approximation 
pi << 1. This transonic approximation was made in equation (18), 
and so forth, only to show the agreement between the common tran- 
sonic limiting value of Ty and T c with the T of reference 30. 

Explicit solution by choice of curvature function. - Before 
choosing a form for the curvature function f = c/c a and thence 

calculating the flow pattern explicitly, a few general remarks, 
which may be helpful in judging the degree of approximation of the 
resulting solution, will be made. The present formulation of the 
airfoil flow problem is in terms of longitudinal distance x and 
local velocity v as independent variables and local streamline 
curvature C and lateral distance y as dependent variables, 
that is, the solution for given stream Mach number M Q is 
required in the form C(x, v), y(x, v) . Various possible pro- 
cedures for obtaining this solution suggest themselves. For example, 
an iterative procedure would be to assume in first approximation a 
curvature function C(x, v), and thence calculate y(x, v) and, 
in particular, the streamlines, using equations (l) to (9). The 
curvature of these streamlines would then yield a second approxi- 
mation C(x, v) and the procedure would be repeated until satis- , 
factory convergence was achieved. Again, a network- type procedure 
would be to set up the continuity integral (l) between n pairs of 
streamlines rather than between the one pair (airfoil and infinity) 
in equation (4). The number n would increase in successive steps 
of the calculation. In each step a curvature function based on the 
curvature values obtained in the preceding step would be used to 
yield by simultaneous solution of the continuity and distance equa- 
tions (1) and (9) an increased number of curvature and corresponding 
distance values for use in the next step. 

The rapidity of convergence of the preceding types of procedure 
could be expected to depend strongly on the initial choice of cur- 
vature function in the first approximation. It is therefore sensible 
for computational reasons to choose an initial curvature function in 
the light of the information supplied by experiment and other theory 
as to the behavior of the streamline curvature in the type of flow 
field under consideration. 
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Only the first approximation is considered herein and free use 
is made of experiment and other theory in setting up an initial 
curvature function. In view of the preceding remarks, this use of 
experiment and other theory should not be construed as inherently 
limiting the degree of rigor in principle attainable in solutions 
by the present approach. This ultimate degree of rigor can be 
expected to be about the same as with existing perturbation or net- 
work formulations of the airfoil flow problem. 

An explicit form for the curvature function c/c a = f(x, z, K), 
equation (13), in the small-perturbation case is now set up. This 
function must satisfy the boundary conditions of the problem; namely, 

at airfoil. 


at infinity. 


V = A 

z = 1 \ 

C =C a 

f(x,l,K) =1 , 


(32a) 


y — » oe> 

V— > o 
z — iO 
C->0 

' f (x,0,K) — > 0 




> 


J 


(32b) 


Furthermore, it may be expected, at least in subsonic flow and at 
the chordwise locations near the maxi-mum velocity or thickness 
regions of the airfoil (which are of most interest), that the stream- 
line curvature will decrease monotonically from the positive value 
C a at the airfoil to zero at infinity. A simple curvature function 
that satisfies these requirements is (reference 29) 

f = z* 


in which the parameter r is a function of x and K to be 
determined. 


( 33 ) 
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Substitution of equation 
tion yields 


Equation (34) gives the local velocity increment A on a symmet- 
ric airfoil in compressible potential flow in terms of the local 
shape parameter YC a , the Mach number parameter p, the specific- 
heat-ratio parameter r c , and the parameter r. 

In order to determine r, it is first noted that to satisfy 
the boundary condition at infinity (equations (32b)), r must be 
positive and the integral for the lateral distance y in equa- 
tion (15) must diverge at z = 0, and therefore r >1. In 
addition, the continuity integral (equation (14)) must converge at 
z=0, which it does for r < 2. Hence r is restricted to the 
range 


(33) into equation (14) and integra- 


2 

jxA 

2-r 


r c A 

3-r 


(34) 


1 < r < 2 


(35) 


In the speed range for continuous potential flow, the flow pattern 
is largely subsonic, as mentioned in the INTEODUCTION . Hence, in 
this speed range the determination of r will be completed (l) by 
assuming r to be independent of K, and (2) by requiring equa- 
tion (34) to yield the (presumed) known velocity distribution 
A i (x) in the incompressible case Mq = 0. Conditions (l) and (2) 

are appropriate to the Prandtl-dauert speed range. In the higher 
speed ranges discussed in subsequent sections in which the flow 
pattern differs very substantially from Prandtl-dauert type, these 
conditions on r will be discarded for new ones. By condition (2) 
and equations (lib) and (lie), equation (34) yields for the deter- 
mination of r 


YC 


a 


_ *i! 

2-r 3-r 


(36) 


The values of r determined by equation (36) turn out to be in the 
range of equation (35) over the maximum- velocity region of the 
airfoil, for which the curvature function (equation (33)) was set 
up. This type of inner consistency will recur often in this report 

The second term on the right side of equation (36) can, in 
principle, be neglected in the small-perturbation case considered. 
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Better results are obtained, however, in applications to airfoils 
of finite thickness ratio (of the order of 0.1) if it is retained. 
Some theoretical explanation for this improvement and hence retention 
of the term in question can be made as follows: The right side of 

equation (34) actually consists of the first two terms of an infinite 
series in the velocity increment A, which is evident for any cur- 
vature function f, from equations (8), (10), and (14). If the 
coefficients multiplying the powers of A are of comparable mag- 
nitude, as in equation (36), it may therefore be expected that for 
small but finite A the first two terms of the series would be a 
better approximation to the complete series them the first term 
alone. 

Thus, finally, in order to calculate a velocity increment A 
at a chordwise location x on a symmetric airfoil at zero lift at 
Mach number M Q by the present method, the shape parameter YC a 

and the low- speed velocity increment A^ at that point are deter- 
mined first. The quadratic equation (36) for r is then solved 
for the value of r in the range of equation (35). (if such a root 
of equation (36) does not exist, the present explicit solution based 
on the curvature function (equation (33)) fails.) While keeping 
this value of r constant as Mq varies and using the values of ji 

and T c corresponding to the desired Mq, equation (34) is solved 

for the desired velocity increment A(x) at free-stream Mach number 
Mq. For convenience in solving the cubic equation (34) for A for 

various values of the other parameters, these other parameters can be 
combined as follows: Let 


A — 



IcA 

d 


X s 


(2-r) 5 Ufoa 
(3-r) 2 li 3 


Equation (34) then becomes 



(37a) 

(37b) 


(38) 


Equation (38) is plotted in figure 3. The most significant feature, 
the double-valuedness of (positive) A for given will be sub- 

sequently discussed in connection with the potential-limit 
phenomenon. 
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Conditions in flow field. - Once the velocity increment A at 
the airfoil 'has been obtained; the variation of local velocity v 
in the field with distance y from the airfoil can he determine 
from equation (15). With the curvature function (equation (33)), 
equation (15) yields 




(r-l)(y-Y)C a 

5 J 


ir- 1 


(39) 


The streamline curvature is expressed in terms of 
tions (33) and (39) as 

-r 
■ r-1 

(r-l)(y-T)C c 

f = Zr « L 1 + 

■'a *- 


y by equa- 


sf 


(40) 


The preceding equations show that both the local velocity and the 
streamline curvature decrease monotonically with y from their 
values at the airfoil to zero at infinity. The actual manner of 
this decrease depends, at a given chordwise location x, on the 
airfoil curvature C a and on the equations of motion through the 

parameters r and A. Actually there is therefore not a great deal 
of arbitrariness in the values of streamline curvature obtained here. 
Once a basic form such as equation (33) is chosen, the actual values 
of streamline curvature in the flow field are determined by the 
equations of motion and in such a way as to satisfy these equations 
on the average or in the large (inasmuch as the integral, rather 
than the differential, form of the continuity equation was used; . 

The insensitivity of airfoil velocity and hence also of the flow 
pattern generally to the choice of a particular form of curvature 
function such as equation (33) is thus plausible. 


A streamline in the flow field can be traced out as follows: 


A streamline is defined by a constant value of 
equation (l), which can be written as 


v in the continuity 

v O 


^o “ 



( 41 ) 
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Eliminating v in equation (41) by the irrotationality relation 
(5) gives the equation of the streamline 

y = S(x,y Q ) (42) 

Thus j Q is in effect the stream function. In the small- 
perturbation case, equation (41) can be transformed as follows: 




7 0 - (y-T) = 


(pv-l)dy 


JTh 


dv 

C 


nA 


= JL 
c Q 


Kj 


ni 


V 


(uu- r v ) 


dv 




z - 


IcA 2 i dz 
z 


.2 

c 


A 

a t- 


(i- z 2 " r ) _ r c A (i-z 3 ~ r ) 

2-r \i 3-r 


Hence, using equation (34), 


2 / 2-r p a z 3 ~ r 

V _ T = z 1£1 _ 

y C a \ 2-r u 3-r, 


(43) 


Elimination of z in equation (43) by use of equation (39) gives 
the equation of the streamline (equation 42). The inclination 0 
of the local stream direction to the x-axis can then be obtained 
as the slope of the streamline. 


It is of interest to note that for r = 3/2, which is in the 
middle of its permissible range (equation (35)), the local velocity 
increment far from the airfoil is given in the Prandtl-Glauert 
range by equations (39) and (16) as 


lim v 
y-> » 


PW 5 


(44) 


This expression agrees with that obtained for a source-sink doublet 
by conventional methods. 
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Illustrative examples. - In order to illustrate tlie preceding 
explicit solution, consider first the Kaplan section (reference 6) 
of thickness ratio t = 0.1. This contour, indicated in figure 4, 
can he expressed in parametric form as (reference 34) 

2Y s t\ = | t (sin cp - | sin 3cp^ (45a) 

2x-l = £ = (l - ^ cos cp + j cos 3Cp (45h) 

where 0 < < P<2ir. When the shape parameter YC a (x), as deter- 
mined hy equation (45), and the incompressible velocity distribu- 
tion Aq(x), as determined by conformal transformation, are sub- 
stituted into equation (36), the parameter r(x) is found to vary 
from r = 1.301 at the midpoint £ = 0 to r = 2 at | = 0.645 
(because A ± = 0 but YC a ^ 0 at g = 0.645). Thus, by 
equation (35), the preceding explicit solution fails outside the 
range 0 < | £|< 0.645. The velocity distributions A(x) at var- 
ious Mach numbers Mq, as obtained from equation (34) within this 

chordwise range, are shown in figure 4. For comparison, the results 
obtained (reference 29) without making the small-perturbation approx- 
imation of equation (10 ) for the density and the results of Kaplan 
(reference 6) are shown. The results of Kaplan were obtained by 
substituting a series for the potential in powers of T into the 
partial differential equation for the potential and successively 
solving the Poisson equations obtained by equating orders of mag- 
nitude in T. The calculation was carried out to the third order, 

3 

T . 


The three sets of distributions in figure 4 are seen not to 
differ greatly except for Mach numbers near the potential limit, 
that is, the maximum Mach number for which the method yields a 
solution. This limiting case will be discussed in the next section. 

As a further illustration, A(Mq) at the maxi mum thickness 
locations of the three T = 0.1 sections - Kaplan, symmetric 
biconvex, and ellipse - as calculated by the small-perturbation 
integral method, the finite-perturbation integral method (refer- 
ence 29), and by various other authors is shown in figure 5. 

The results of Hantzsche (reference .7) for the ellipse were 
obtained by a method similar to that used by Kaplan. The calculations 


12^5 
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were carried to the third order in a thickness-ratio parameter, which 
in successive orders equalled T, T 2 , and T 3 + 0(t^ log t), 
respectively. 

In the method of Gothert- and Kawalki (reference 17), the influ- 
ence of compressibility was represented by a' finite number (12 per 
quadrant) of sources and sinks (Bayleigh- Janzen idea) located in the 
incompressible-flow field, as corrected by the Prandtl-Glauert rule, 
outside a given airfoil. The source-sink strength used for calcu- 
lating velocity increments caused by compressibility was computed 
by iteration from the flow-field velocities themselves so as to 
satisfy continuity. As many successive steps in the calculations 
( ~ 10 near the potential limit) were made as required for obtaining 
a self-consistent compressible-flow field that could be regarded as 
an incompressible-flow field, with the given spatial arrangement of 
sources and sinks superimposed. 

Finally, the velocity increments given by the Prandtl-Glauert 
and Karman-Tsien rules, (reference 35) are included in figure 5. 

It can be seen from figure 5 that the results of the present method 
differ from those of the other methods by about the same amounts as 
the results of the other methods differ amongst themselves. 

The Prandtl-Glauert rule has an interesting connection with 
the present method. This rule results when the second term in 
equation (34) is negligible relative to the first term, as Is the 
case in the subsonic speed range for small perturbations. The tran- 
sonic term is also negligible, however, for higher speeds or larger 
perturbations if r ~ 2. As previously noted, r — >2 at those points 
on the airfoil where ~ 0 but YC a f 0. Near these (locally 

subsonic) locations the Prandtl-Glauert rule may therefore be expected 
to apply. On the other hand, the maximum velocity increments on air- 
foils of the type illustrated in figure 5 turn out to vary with Mach 
number in the subsonic range in a manner approximating the Karman- - 
Tsien rule. (See also reference 29). Inasmuch as the KArm&n-Tsien 
rule itself approaches the Prandtl-Glauert rule for small velocity 
increments, the present method therefore gives results in the sub- 
sonic speed range not greatly different from those obtained by the 
Karman-Tsien . rule . 

Potential-limit phenomenon. - In computing a velocity distri- 
bution on an airfoil at a fixed free- stream Mach number Mq by 
equation (34) as plotted in figure 3, a portion of the airfoil con- 
tour is put into correspondence with a portion of the lower branch 
of figure 3. For example, on the Kaplan section at Mq = 0, the 
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region of airfoil 0.178 < x <0.823 corresponds to the portion AD 
in figure 3, such that traversing the stretch 0.178 < x <0.5 on 
the airfoil corresponds to traversing the stretch AD and traversing 
the stretch 0.5 < x < 0.823 on the airfoil corresponds to trav- 
ersing the stretch DA. 

As the Mach number increases, the same region of the airfoil 
corresponds to a greater portion of the lower branch of figure 3, 
the point A corresponding to x = 0.178 and 0.823 remaining at the 
origin but the point D corresponding to x = 0.5 moving away from 
the origin. For example, the fixed airfoil location x = 0.5 corre- 
sponds at Mq = 0.5 to point D^_, at Mq = 0.75 to point and 

at M Q st Mq j = 0.807 to the turning point B between lower and 

upper branches in figure 3. Thus, at a certain lowest free- stream 
Mach number Mq^ a point X 2 on an airfoil comes into corre- 
spondence with the turning point B in figure 3. For higher but 
still subsonic free-stream Mach numbers, the stretch AB remains 
fixed at its maximum extent but now corresponds to a di mini shing 
region of airfoil, such that the gap on the airfoil for which no 
solutions for the velocity are yielded by equation (34) or figure 3 
spreads out from x^ as M Q increases in the range M q ^ < M q <1. 

For example, on the Kaplan section the forbidden region at Mq = 0.85 
is 0.35 < x <0.65. Thus, Mq^ is the maximum free-stream Mach 

number for which continuous potential -flow solutions are obtainable 
by the present method. Herein is termed "the potential-limit 

Mach number," and the corresponding flow pattern, "the potential- 
limit solution." 


Quantitatively, the turning point B corresponds to the condition 
dx/dA =0 in equation (38) or, equivalently, the condition 
(du/dA) rj ,YC a = const = 0 equation (34). Hence, at the potential 

limit (denoted by the subscript l). 


\ 



(46a) 


or, by equation (37) 




2 

3 


5-r 

2-r 


(46b) 


l 


(47a) 
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_ 21 (2-r) 3 YC 
" 4 (3-r) 2 S 


(47b) 


The procedure for obtaining the potential- limit Mach number for a 
symmetric airfoil at zero lift is then to plot the right side of 
equation (47b) as a function of chordwise location x. The potential- 
limit chordwise location x^ is that where the right side attains a 

maximum value. This maximum value yields, by equation (47b), a value 

3 3 

of p. /r that, by equations (lib) and (lie), determines the 
potential-limit Mach number M^ • The corresponding velocity- 
increment Aj is obtained from equation (47a). The velocity dis- 
tribution over the rest of the airfoil at Mach number M is 
determined as usual from equation (34) or (38) using figure 3. 

The potential-limit velocity distributions for the Kaplan sec- 
tion by the small- and finite-perturbation integral methods are 
shown in figure 4. The potential-limit Mach numbers by these two 
methods are 0.807 and 0.843, respectively. Kaplan (reference 6) 
estimates from the first three terms of his expansion for the poten- 
tial that the entire series would diverge beyond Mq = 0.83. Gothert 

and Kawalki (reference 17) also locate a limiting Mach number above 
which their method does not converge to a self-consistent solution 
regardless of the number of iterations performed. Although the var- 
ious determinations of limiting Mach number M and the corre- 
sponding velocity increment , indicated by the circled points 

in figure 5, do not agree extremely closely, nevertheless, consider- 
ing the diversity of the methods used, the inference appears war- 
ranted that a unique, upper- limiting, subsonic free-stream Mach num- 
ber exists beyond which a continuous potential flow past a given 
airfoil cannot be derived. 

Various properties of the potential-limit flow pattern can be 
derived by the present method, which not only make it easy to under- 
stand why this limiting flow pattern should exist but also connect 
it with the so-called limiting-line phenomenon encountered in the 
•hodograph method. (See, for example, reference 36.) In the first 
place, at the potential limit the local Mach number Mj at the 
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chordwise location x^ on the airfoil must he greater t han 1, 
hence a local supersonic region must exist in the flow field. By 
equations (27) and (47a), 


M, 



» l 


4 (3-t\ 

3 r c \2-tJ 


- 1 


(48) 


In the permissible range (35) for r, the 'right side of equa- 
tion (48) is positive down to at least Mq j =0.63. For the thin 

airfoils contemplated here, Mq j > 0.63 and hence > 1, as 
indicated by the constant local Mach number contours in figure 5. 

Secondly, the velocity gradient is finite and suffers a 

finite discontinuity 1 at x = x^ when M Q = Mq l . Differentiation 
(denoted by a subscript) of equation (34) with respect to x yields 



r c A 5 

( 3 -r ) 2 


( 49 ) 


The quantity in parentheses on the left side of equation (49) is 
zero for M Q = M Q ^ , by equation (47a) . The right side of equa- 
tion (49) is zero at x = x ? for Mq = Mq as can be seen from 
equation (47a) and from the potential-limit condition that 
in equation (47b) be a maximum at x = x^; that is, that its 
derivative by x is zero at x^ . Thus, x in equation (49) is 
given by an indeterminate form at Mq = Mq^, x = x ? . Application 

of 1’ Hospital’s rule to this indeterminate form and use of equa- 
tion (47) yields a quadratic equation for A x . For the usual air- 
foil shapes, the two roots of this quadratic are real and of oppo- 
site sign. (For two-way symmetric airfoils such as in fig. 5, 
the two roots are numerically equal.) Hence A undergoes, in 

general, a finite discontinuity at x = x^ for Mq = Mq ^ of an 

amount equal to the difference between these two roots. This dis- 
continuity in velocity gradient is indicated for the Kaplan section 
in figure 4 at x = 0.5 and = 0.807, M Q 1 = 0.843 by the 

small- and finite-perturbation methods, respectively. If 
(YCa^xx ~ T x = r xx = 0 at x Mq = M^, then also A x = 0 

and no discontinuity occurs in A x at the potential-limit point 
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in this case. The preceding statements for also hold for the 

fluid acceleration on the airfoil inasmuch as the acceleration is 
given by (l+A)^. 

Third, hy eliminating y between equations (39) and (43) and 
differentiating the result by x, an expression for the veleoity 
gradient (V 2 ) s along a streamline in the flow field is obtained 

that is a linear function of A^, (YC a ) , r x , and C a x , with 

coefficients dependent on A, YC & , r, and C & . Furthermore, the 

slope of a constant velocity, or constant local Mach number, line 

(y x ) , multiplied by the local streamline curvature C, equals 
M 

(v x ) . Hence (V x ) , the local fluid acceleration along a stream- 
s s 

line, and (y*-) w are associated with A at a given x in being 
finite or infinite together. 


Fourth, the present method yields a second flow pattern at and 
above the potential limit, which overlaps the first. This pattern 
is given by the upper branch of the cubic curve, figure 3; for, at 
the potential limit, the turning point B is Just reached from along 
the lower branch AB upon traversing the airfoil to its potential- 
limit chordwise location x, . In retracing the airfoil in the oppo- 
site direction, the upper branch BC could be used. For Mq > Mq^ y 
instead of a single chordwise location x^, two chordwise locations 


Xj ,i and x -^2 correspond to the turning point B, as previously 
mentioned. Between the two lines x = x ^i arL d. x = x^g “thus 


a region of space for which no solution can be found by the present 
method. Outside this region two overlapping flow patterns are 
derivable, corresponding to the two branches in figure 3. At the 
lines x = x~i x = x^ 2 > the velocity gradient is, in gen- 

eral, infinite. (The left factor in equation (49) is zero by equa- 
tion (47a), but the right side is not zero, because 2 ad- 


equation (47b) is here not a maximum.) These results and those of 
the preceding paragraph show that the lines x = x^-j_‘ and x = x^^j 

which are potential lines in the present approximation, are envelopes 
of constant-velocity, or constant Mach number, lines. Furthermore, 
it is plausible that the lines x = x^]_ and x = x^2 are also 

envelopes of Mach lines or characteristics; a point of infinite 
velocity gradient in a supersonic region requires a spacing of Mach 
lines closer by an order of magnitude than does a point of finite 
velocity gradient, and an envelope of Mach lines satisfies this 
condition. 
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The preceding properties indicate that z = Zj i am ^ 
z = Zj ^2 “a? he limiting lines in the sense of the hodograph method 

and that the potential-limit solution is the boundary solution between 
continuous potential-flow patterns, for Mq < and flow patterns 

with limiting lines bounding forbidden regions, for Mq > 

inference is supported by a closer examination of the limiting-line 
phenomenon as follows: 


A limiting line in a flow field is defined as the locus of 
points at which the Jacobian J = 5($,''lO/d(v,0) from the hodograph 
(v,0) -plane to the physical (0 , ¥) -plane vanishes. (Constant $ 
and ¥ define the potential, and streamlines, respectively, in the 
physical plane.) Consider the reciprocal of this Jacobian 


_-l 9(v,0) 9v 90 c>v 90 


(50) 


The differential equations of motion in intrinsic form are (refer- 
ence 35) 


90 _ (M 2 -!) 9v 
- pv 


(51a) 


96 p 9v 
= v^¥ 


(51b) 


Substituting equations (51) into equation (50) and noting that 
d$ = vds and dUT as pvdn yields 


j- 1 . pv3 ** = (if-i) (gf - (||) 2 (*» 

where 9v/9s is the velocity gradient along a streamline and -90/9s 
is the streamline curvature. 

The variation of j"-*- along the t = 0.1 biconvex airfoil, 
determined by equation (52) and the velocity distributions previ- 
ously calculated, is shown in figure 6 for various values of Mq. 

Figure 6 shows that the present method yields finite values of 

in the range M Q < M The only discontinuity in behavior 
occurs at z = z^ =0.5 at which point undergoes a finite 
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discontinuity for Mq = Mq ^ (“because, as previously discussed, 
Sv/Ss undergoes here a finite discontinuity as a function of Mq). 

It may be inferred, in general, that < 0 everywhere for 
Mq < Mq i. For, equation (52) can be solved for (Bv/d0) s as 

(subscript s denotes values along a streamline) 



For = 0, equation (53) yields the Prandtl -Meyer rate of change 

for (Sv/d0) s in supersonic regions. For equation (53) 

yields a numerically greater rate. How, in a local supersonic region t 
embedded in a subsonic flow field, it is known that (Sv/d0) B must 

be numerically less than the Prandtl -Meyer value (references 37 and 
30) because the sonic boundary reflects Mach waves of opposite sign 
to those incident upon it. Flows having finite fluid accelerations 
and streamline curvatures and containing local supersonic regions 
must therefore have a finite < 0 everywhere (in subsonic 

regions j”^-< 0 by equation (52)). Hence, also j < 0 within 
such flow fields. The preceding discussion may be regarded as an 
indication of a simple proof of a theorem of Friedrichs (reference 39) 
to the effect that if a continuous potential-flow solution v(x,y,M Q ) , 
0(x,y,MQ) varies continuously with Mq and approaches a limit as 
Mq approaches some limiting value, then this limiting solution 
itself cannot contain any points at which j = 0; that is, a 
limiting line. 

The relatively small region of positive values of j - - 1 - indicated 
in figure 6 for Mq < Mq^j must therefore be regarded as due to, and 

indicating the extent of, the approximations in the present method. 

A plot of j-1 would appear to be a useful indication for assessing 
and possibly improving the accuracy of transonic velocity distribu- 
tions. The finite discontinuity in J - 1 (Mq) at x = x^ for 

M 0 = Mq j may well exist (in theory). The present method thus 
ill umina tes the theorem of Friedrichs. 

For Mq> Mq j, figure 6 shows that j - -'- — > co, hence j— >0 
at a point x^ ]_(Mq) and at the symmetrically disposed point 
x l 2 = 1-x 2,l (because here dv/ds— >», as previously noted). 


t 
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Thus, formal continuation of the present method to M Q > M yields 
discontinuous solutions containing forbidden regions x^ ^ < x < x^ ^ 
hounded by limiting lines, as previously inferred. 

The reason for the existence of two values of X for a given X 
in figure 3, and hence two overlapping flow patterns, of which one 
corresponds to either subsonic or slightly supersonic conditions 
and the other to highly supersonic conditions, is essentially sim- 
ilar to the reason for the existence of two possible velocities for 
a given flow area in one- dimensional flow, as in a converging- 
diverging nozzle. Indeed, the present method, which makes use of 
the equation of continuity in quasi- one- dimensional or hydraulic 
fashion, shows rather directly that the potential-limit solution 
occurs for essentially the same reason that a converging-diverging 
nozzle chokes when a sufficiently high but subsonic Mach number 
occurs ahead of the throat (reference 29) . (The mass-flow density 
pv passes through a maximum at a local Mach number of 1.) The 
manner by which steady-flow patterns can nevertheless exist for 
the subsonic free- stream Mach numbers between M Q ^ and unity, in 

spite of the analog of this choking phenomenon (which prevents higher 
subsonic Mach numbers from being obtained ahead of the throat in 
nozzles) is considered in the next section. In particular, it 
will be seen that the purely mathematical extrapolation made 
herein leading to the limiting lines x = } i an< ^ x « x^^j 

has no physical significance in connection with the shock that 
actually appears at these Mach numbers. 


LOWER TRANS OHIO ASYMMETRIC FLOW WITH 
TERMINAL SHOCK 

Velocity distributions. - In this section, asymmetric flow 
patterns with a single terminal shock will be derived for free- 
stream Mach numbers from the potential- limit Mach number to unity. 

The analysis, particularly that connected with the terminal shock, 
is more uncertain in this range than in the neighboring speed ranges 
because the terminal shock is no longer at most a boundary of the 
region of interest, but exists within it. 

The procedure for obtaining a symmetric solution at the 
potential-limit Mach number, derived in the preceding section, indi- 
cates how asymmetric solutions may be obtained at higher Mach numbers. 
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Thus, in determining the airfoil velocity increments at the 
potential-limit free- stream Mach number as a function of chord-wise 
location, the lower branch of the curve in figure 3 was traversed 
from A to B, the potential-limit point B being reached, say, at the 
midchord location for the symmetric biconvex: airfoil. In going from 
midchord toward the trailing edge on the airfoil, the curve was 
traced back from B to A, giving the symmetric solution. Because the 
turning point B at midchord has just been attained, however, it is 
possible to continue around and have part of the upper branch BC 
correspond to the rear part of the airfoil. The chordwise location 
x = Xj in the present speed range therefore has the property of a 

branch point, as previously explained. This property maintains the 
analogy with conditions at the throat of a converging- diverging 
nozzle . 

The preceding procedure yields an asymmetric solution at the 
potential-limit Mach number. In order to obtain a solution, and in 
particular an asymmetric solution, without forbidden regions at 
higher subsonic Mach numbers, the basic cubic equation (34) must be 
modified. The modification made herein is to introduce a nonvan- 
ishing, lateral streamline displacement 6 at infinity (fig. 2). 

Such a displacement, of the order of the airfoil thickness, is 
clearly necessary for steady-flow patterns to exist near sonic speed 
if the free- stream flow is to "get by 1 ' the airfoil In accordance 
with steady- state continuity relations. (That a steady- flow solution 
is yielded by the equations of motion at precisely sonic speed has 
been demonstrated by Guderley, references 38, 40 and 41.) The stream- 
line at infinity can still be considered straight and parallel to the 
x-axis, relative to the airfoil slope, if its lateral displacement 
5 is assumed to occur gradually over an infinite distance along the 
streamline. Thus, 6 will not be a function of airfoil chordwise 
location but will be constant for the entire airfoil. It will, how- 
ever, be a function of free- stream Mach number. The introduction of 
a nonvanishing, lateral streamline displacement 5 in the present 
speed range M Q < M Q ^ < 1 does not conflict with the proof pre- 
viously indicated that 6 vanishes for subsonic-flow conditions and 
no wakes at infinity; the simultaneous introduction of a shock' ter- 
minating an asymmetric local supersonic region implies a shock wake 
of slightly decreased total head extending to infinity. Such a wake 
provides a mechanism for producing a nonvanishing displacement S 
similar to that provided by an ordinary subsonic wake produced by 
boundary- layer action. 
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By subtracting a suitable displacement function 6(M q ) from 

Y in the basic cubic equation (34), both branches in figure 3 can 
now be correlated continuously with the airfoil contour in the 
speed range M < Mq < 1 to yield an asymmetric solution. In 

order actually to calculate this solution, the variation of the 
curvature parameter r with x and K, or its equivalent, must 
first be determined. To this end, the assumption is first made 
that the potential-limit condition of equation (47), with Y 
replaced by Y -5, will apply throughout the speed range 
Mq 2 < Mq < 1 to the same chordwise location x = x^ first found 

to require the potential-l im it condition. In more geometric terms, 
if the basic cubic equation (34) is viewed as a surface in the A , 
x, and Mq space, then the potential-limit condition, which has 

been regarded as determining an absolute maximum of Mq on the 
projections of this surface in the A, M q plane, can therefore 

be regarded as determining the saddle point of this surface. Hence, 
in the A, x plane, the potential-limit point is a branch point. 
Secondly, it is noted that the velocities obtained along the rear 
portion of the airfoil, corresponding to the upper branch in fig- 
ure 3, and using the subsonic values of r previously determined, 
become much greater than those corresponding to a Prandtl-Meyer 
distribution. As pointed out by Tsien and Fejer in reference 37, 
however, the velocity in a local supersonic region should be less 
than that in a Prandtl-Meyer expansion through the same angle of 
turn, the difference being due to compression waves reflected from 
the sonic boundary and striking the airfoil. Inasmuch as these 
compression waves will probably be of small influence compared to 
the expans ion waves emanating from the airfoil for the size and the 
shape of local supersonic regions to be expected in the present 
speed range, it will be assumed that the supersonic portion of the 
velocity distributions to be derived will be Prandtl-Meyer dis- 
tributions. (it would not be difficult to obtain more accurate 
velocity distributions in local supersonic regions on the basis of 
the shape and extent of the supersonic region determined by the 
first approximation considered herein. ) The Prandtl-Meyer velocity- 
distribution formula can be written in the transonic small- 
perturbation limiting case as (reference 30) 

_ 2/3 

m 2 -i = u + 2rA= 3r(e-e s ) ( 54 ) 


vhere 


p = -p = m q 2 -i 


(55) 
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and 0-0 e is the (positive) angular change, in. radians, of slope of 
airfoil contour from the value 9 e at which the local Mach number 

is 1 to the value 9 at which it is M. Throughout this section 
and the nest, for simplicity in the calculations, all T's that 
arise, such as the of equation (28), which should appear in 

equation (54), are replaced by F = 1.2 unless otherwise mentioned. 
Hence, the calculations are limiting calculations for Mq = Ij that 
is, the transonic approximation is made more strictly than in the 
preceding section, where the functions F c (Mq) and were 

used. 


By assuming a Prandtl-Meyer distribution, the problem of cal- 
culating the supersonic portion of velocity distributions in a given 
speed range is reduced to determining one point of this distribution 
(that is, velocity and airfoil location) at each free-stream Mach 
number Mq in the range, which is necessary to evaluate the constant 

0 e in equation (54). This one point will be chosen as the velocity 
at the potential- limit, or branch-point, location x = x^ in the 

present speed range. . It is convenient to carry out the determination 
with the basic cubic equation (34), as modified by inclusion of 5, 
and the potential-limit condition of equation (47) in explicity 
similarity form, thus: 


where 


(Y-5)C a ^ kl 2 l 5 
-.2 ” 2-r “ 3-r 


^ -s 




(Y-B)C. 


3 \2-r J 


% 


(3-r) 

,3 


27 \2-r/ ^ 3 ' r ) 


(56) 

(57a) 


(57b) 


K = n/(rr ) 2/3 

L = TA/irT) 2 ^ 


(58a) 

(58b) 


Explicit calculations will be made for the thin symmetric biconvex 
airfoil. This airfoil contour is given, for small thickness ratio 
T, by 
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| s 2x-l 


9 

2T 


(59a) 


T (1-| 2 ) 


(59b) 


The left side of equation (56) becomes, by equations (59) 


(*-6)C a 

2 



(60) 


The potential- limit location for the symmetric biconvex airfoil 
is Xj = 0.5 or £2 = 0* Inasmuch as a unique steady- flow pattern 

at Mq = 1 is assumed, which is continuous with the flow patterns 
existing on either side of Mq = 1, the velocity at £ = = 0 

in particular must have this property. In the next section, the 
1/3 

value I^ = 2 ' is derived for the velocity parameter at 
|=|=0 at K = 0 or M 0 = 1, by an analysis for M Q > 1. 

In order for equation (57a) to yield this value also, it is nec- 
essary that r-4 2 as K — >0 and in such a way that 


lim fe ■ J73 = 1 - 890 < 61 > 

' K — >0 2 ' 

At sonic speed, a nonvanishing streamline displacement at infinity 
8 is thereby obtained, given by equations (57b), (61), and (60) 
with £ = 0 as 


28/t = 0.5 (62) 

This value of 5, half the maximum airfoil ordinate, is of the 
right order of magnitude inasmuch as at precisely sonic speed a 
streamline displacement at infinity somewhat greater than the max- 
imum airfoil ordinate is clearly required by steady- state continuity 
considerations . 

The velocity at the potential-limit chordwise location £ 2 = 0 
can now be obtained as a function of M Q by an interpolation pro- 
cedure. The value of K/(2-r) at Mq = = 0.831 and 

£=£ 2=0 is 1.60 from the results of the preceding section. This 
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value does not differ greatly from that at Mq = 1 given by equa- 
tion (61). Accordingly,- a linear variation of E/(2-r) at £ = 0 
as a function of K can be assumed in this range that yields r(K) 
and thence Iq (K) and 8(E) at | = 0, by equations (57) and 

(60). The values Lj (E) at | = 0 enable Prandtl-Meyer velocity 

distributions to be determined over the airfoil at each Mach number 
M Q corresponding to K by equation (54). 

As a possible refinement, the velocity distribution over the 
front portion of the airfoil can be interpolated between the Prandtl- 
Meyer distribution at Mq = 1 and the potential-limit distribution 

at M 0 = Mq j found in the preceding section. The interpolation is 

made at each chordwise location x between the values of E/(2-r) 
given by equation (56) using equation (62) and the values of L 
given by the Prandtl-Meyer distribution at Mq = 1, and the values 

of K/(2-r) at M . = 0.831 determined in the preceding section. 

o, L 

The resulting values of r(x,E) and 5(E) determined from equa- 
tion (57b) then yield intermediate velocity distributions by equa- 
equation (56). 

The asymmetric-type velocity distributions determined as Just 
outlined are so far incomplete in two basic respects. First, the 
subsonic portion from the leading edge to the point where free- 
stream velocity is first reached locally is lacking. This lack, 
as explained in the preceding section, is mainly due to the inade- 
quacy of the simple curvature function (equation (33)) in this 
region. For the purpose of calculating pressure drag in the tran- 
sonic range, this portion is not needed. (See section PRESSURE DRAG; 
USE OF MOMENTUM INTEGRAL.) The second and more important omission, 
the terminal shock, is discussed next. 

Terminal shock. - The terminal shock is observed to form the 
downstream boundary of the local supersonic region and to move 
toward the trailing edge as the free- stream Mach number increases 
toward unity. Assume for the present that this shock is normal 
rather than oblique. This assumption appears valid, at least under 
the experimental conditions reported in reference 2, in which the 
measured pressure rises across the shock and outside the region of 
shock-boundary-layer interaction corresponded fairly closely to 
those of ideal (that is, Rahkine - Hugoniot ) normal shocks. 

The velocities across an ideal normal shock are related by 
Prandtl’s equation, derivable from the Rankine- Hugoniot relations. 
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V 1 V 2 = a * 2 ( 63 ) 

where a* is the speed of sound at sonic speed M = 1 in an 
isentropic flow. In the small-perturbation transonic case, equa- 
tion '(63) becomes, for the shock element near the airfoil 
(appendix C), 


+ TA 2 = p (64) 

In terms of the local Mach numbers = aJ 1 + p-^ and Mg = /\J 1 + pg 

upstream and downstream of the shock, respectively, equation (64) 
becomes by equations (27) and (28b) 

ill + p 2 = 0 (65) 

Equations (64) and (65) are unsatisfactory for determining terminal- 
shock locations; because first if no boundary layer is assumed and 
subsonic values of Ag are determined either by a subsonic rule as 

in reference 34, or equivalently by the subsonic results of the pre- 
ceding section, then the chordwise locations at which such subsonic 
values of Ag will "match" the supersonic values of A^ previously 

determined to satisfy equation (64) or (65) are within about 20 per- 
cent of the chord length of the trailing edge for the symmetric 
biconvex airfoil. This result would probably not be substantially 
altered by a more exact calculation of viscosity- free flow patterns 
with single terminal shock, such as by the relaxation method of, 
reference 18. Experiment shows, however, (references 3 and ,4) that 
a well-developed shock exists forward of the aforementioned location. 
On the other hand, if a boundary layer is assumed that, in traversing 
the shock, separates so completely as to make the pressure downstream 
of the shock approximately equal to free- stream pressure, then TAg 

(the velocity outside the boundary layer) drops out of equation (64). 
The shock locations determined by TA^ = p, however, move forward 

instead of rearward as sonic free- stream speed is approached from 
below. An intermediate situation, in which Ag could be determined 

from velocity distributions calculated to allow for the presence of a 
boundary layer that has interacted with the shock, might yield sat- 
isfactory terminal-shock locations by equation (64). Such a calcula- 
tion of Ag has not yet been attempted because of” the difficulties 

of calculating the shock -boundary- layer interaction and Its effect 
on the downstream velocity distribution. 
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Some understanding of, the rearward movement of the terminal 
shock can be had,- however, from a modification of equation (64). 
Assume that the actual pressure rise across the normal shock differs 
from the Bankine-Eugoniot, or ideal value, by an amount Ap g . This 

difference can be caused inside the region of shock -boundary- layer 
interaction by friction processes, in which case the actual pressure 
rise is less than the ideal value. A pressure drop Ap_ can also 

D 

occur outside the boundary layer, according to references 2, 18, 
and 19. This drop occurs over a short distance just downstream of 
the (ideal) normal shock and exists because of the necessity for 
the flow to make a quick reversal of curvature from the negative 
value given it by the shock to the positive value required by the 
airfoil contour. Farther out along the shock, this pressure drop 
becomes a pressure rise for the same reason (reference 2). In any 
event, this type of pressure change must physically "join on" con- 
tinuously with that in the boundary- layer region. In view of the 
preceding discussion, a shock pressure-drop coefficient c s , rep- 
resenting the over-all departure of a normal-shock element from 
Bankine-Hugoniot behavior, can be defined in the small-perturbation 
transonic case as 


c s = 


Ap e (P 2 -P 0 ) 


Pi v i 


ide al 


(P2'P 0 ) 


actual 


2 p o V o 


2 p o V o 




( 66 ) 


If the element of shock front under consideration is near the air- 
foil, as is the case here, then Pg in equation (66) is the pres- 
sure coefficient at the airfoil just downstream of the' shock and 
A 2 is the Bankine-Eugoniot velocity that would exist just down- 
stream of the shock if c s were zero. Substitution for Ag from 
equation (64) yields 

- 2 (jl - f) - P 2 ( 6? ) 

The additional assumption is now made that the pressure on the 
downstream side of the shock equals free-stream pressure, or 
Pg = 0. (in a conversation with the author, von Karman mentioned 

having used this condition some years ago.) This condition, in 
general, requires separation of the boundary layer as it traverses 
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tlie shock. If c 6 were known as a function of M-j_ (and of a 
suitable Beynolds number at the shock if c s represents a fric- 
tional effect), shock locations could be determined from equa- 
tion (67) with P£ = 0. In the absence of such knowledge, however, 

an assumed rearward shock movement can be used to calculate c B . 

Thus assume the following boundary conditions on the terminal- shock 
location: (a) At the potential-limit Mach number the ter- 

minal shock occurs at the potential-limit chordwise location x^ 

(x^ = 0.5 for the biconvex airfoil); and (b) At u a free-stream Mach 
number of unity, the terminal shock is at the trailing edge, x = 1. 

These boundary conditions are, of course, approximate. The 
terminal shock appears to be formed by the ever-present random pres- 
sure disturbances that are continually propagating through the flow 
field at locally sonic speed. Some of these disturbances, traveling 
upstream (reference 1), encounter the locally supersonic region and 
tend to pile up; that is, the terminal shock appears to form as an 
envelope of characteristics in the x,y,t space (t is time) rather 
than as an envelope of characteristics in the x,y plane (limiting 
line). It may be necessary, however, for the symmetric-type locally 
supersonic region to reach a finite size before the unsteady dis- 
turbances can coalesce to form a steady terminal shock (reference 4). 
The boundary condition (a) idealizes this process by assuming no 
shock and a symmetric solution up to the potential-limit Mach num- 
ber, thence a sudden change-over to an asymmetric solution with 
steady terminal shock. Presumably, the disturbances that form the 
terminal shock have a tendency to pile up as far upstream as possible. 
The terminal shock could hardly form upstream of the potential- limit 
location, however, because of the "choked" condition at this location. 
(See preceding section. In this connection compare the experimental 
data and the semiempirlcal analysis of reference 24.) As regards 
boundary condition (b), there axe no experimental data on terminal- 
shock location and structure at Mach number unity. In the super- 
sonic speed range, this shock becomes the oblique shock, which does 
start close to the trailing edge. The schlieren photographs of 
Liepmann (reference 3) suggest that the terminal shock will reach 
the vicinity of the trailing edge at sonic speed but will also be 
curved in the direction ultimately required in the completely 
supersonic regime. 

Boundary condition (a) and equation (67), with P 2 = 0, now 
yield for the T = 0.1 biconvex airfoil, at Mq = M 0 ^ = 0.831 
and A 1 =s Aj = 0.389 (values shown in fig. 5), the shock pressure- 
drop coefficient c s = 0.26. At Mq = 1 and x = 1, A^ = 0.47 
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(from the following section and fig. 7) giving c s = 0.94. The 
assumption of a linear variation of A^ -with (i between the pre- 
ceding end values (the extreme values of do not differ greatly) 

determines, from the velocity distributions previously calculated, 
the variation of terminal- shock location x s with Mach number Mg 

(fig. 8). Typical velocity distributions as thus completed by 
location of the terminal shock are shown in figure 7. (The velocity 
increment A is related to the pressure coefficient P by 
P = -2A. By the assumption P 2 = 0, A is therefore indicated as 

zero behind the shock. This A refers to conditions outside the 
boundary layer.) 

Other shock conditions have been advanced. Tsien and Fejer 
(reference 37) regard the shock as an oblique shock producing max- . 
imum outward deflection of the streamline outside the boundary 
layer. Theodorsen (reference 42) regards the shock as the oblique 
shock that produces locally sonic speed on its downstream side. 

In the small-perturbation transonic case, these shock conditions 


lead to equations similar in form to equation (64): namely 

(appendix C), 


Maximum deflection: 

rAj/3 + ta 2 = 2/3 ji 

(68) 

Crocco’s point: 

r Aj/7 +TA 2 = 4/7 n 

(69) 

Sonic point: 

r a 2 = 1/2 n 

(70) 


An additional shock condition has been included in equation (69). 

This condition is the oblique shock for which the streamline cur- 
vature on the downstream side is zero (if the oblique shock has 
finite curvature, reference 43). It is defined by Crocco’s point 
on the shock polar (appendix C). 

Prom the similarity between the preceding equations and equa- 
tion (64), essentially the same objections as regards ability to 
predict shock location that have been made to equation (64) can be 
made to the preceding equations. In addition, as has been mentioned,, 
measurements indicate the terminal shock to be normal under the 
experimental conditions of reference 2. Making the same mod- 
ifications, however, in equations (68) and (69) as were made in 
equation (64), that is, introducing a shock pressure-drop coefficient 
c s and assuming P 2 = 0, yields, with the same shock boundary con- 
ditions as before, practically the same variation of shock location 
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with Mach number Mq but different ranges of variation of c s . 

These conditions are listed in table I. (Equation (70) does not 
yield a solution for shock location inasmuch as, unlike the previous 
equations, it does not relate conditions across the shock but speci- 
fies only the velocity on its downstream side. Equation (70) there- 
fore depends in an essential way on the knowledge of the velocity 
distribution downstream of the shock.) 


TABLE I 


Terminal- shock condition 

Shock-pressure-drop 
coefficient, c s 

EH 

M c = l 

Normal 

0.26 

0.94 

Maximum deflection 

- .09 

.31 

Crocco’s point 

- .18 

.14 

Sonic point 

- .26 

0 


The pressure measurements of reference 2 indicate values of 
c s of about 0.1 just outside the boundary layer and of about 0.3 

at the surface of the airfoil at an effective free- stream Mach number 
probably much closer to the potential- limit Mach number than to 
unity. The surface-pressure and wake survey measurements of refer- 
ence 44 can be interpreted as indicating values of c 8 up to about 

0.5 in the boundary- layer region. These estimates of c s would 
appear to favor the normal-shock condition (table I), particularly 
inasmuch as the negative values calculated in table I are probably 
unrealistic. The actual shock-boundary- layer structure near the 
airfoil is, of course, considerably more complicated than that of 
a single straight shock, whether normal or oblique. The introduction 
of a coefficient such as c s to represent an over-all departure 
from ideal single-shock behavior, however, may possibly be useful 
•practically. 

The reason for the rearward movement of the terminal shock with 
increase of Mach number is not necessarily indicated by the pre- 
ceding correlation between terminal- shock location and pressure-drop 
coefficient c s . It would seem more likely that a given terminal- 
shock location resulting from other, more powerful, dynamical 
requirements should cause a particular value of c s to exist, 
rather than conversely. These other dynamical requirements may be 
connected with the necessity for the streamline displacement at 
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infinity S(M q ) to increase with Mach number in the present speed 

range. An increasing 8 could he produced hy a more extensive 
and intensive shock -boundary- layer wake, which in turn would require 
a longer and stronger terminal shock, such as would occur at the 
higher local Mach numbers farther to the rear of the airfoil. 

It need hardly be emphasized that the calculations of this 
report, particularly those in this section, are of a provisional 
and exploratory character. At the present stage, not much more has 
been rigorously demonstrated than that the methods used herein can 
yield results of the right order of magnitude. 

Finally, it is of interest to note that the condition used in 
this section, that shock back pressure shall equal free-stream pres- 
sure, can be regarded as a special case of a more general condition; 
namely, that the shock back pressure shall be as close as it can get 
to free-stream pressure. This more general condition is discussed 
in appendix D. • . 


UPPER TRANSONIC FLOW WITH DETACHED SHOCK 


Analysis. - The upper transonic or detached shock regime extends 
from free-stream Mach number 1 to that supersonic Mach number Mq^ 

at which the shock ahead of the airfoil becomes attached to the 
leading edge, assumed sharp- edged. In the transonic small- 
perturbation case, Mq a is given by (reference 30 and appe n di x C), 


M 2 -l 
o,a 


(OrD 


273 2 2 73 


= 1.890 


(71) 


where 0^ is the semivertex angle at the leading edge (fig. 9). 

Calculation will now be made of the movement of the head shock 
(fig. 1) and of the supersonic part of the airfoil velocity dis- 
tribution as a function of Mq. The terminal shock is assumed to 

start from the trailing edge (see preceding section) so that it need 
not be considered in the analysis. As in the preceding section, the 
supersonic part of the velocity distribution on the airfoil is 
assumed to be a Prandtl-Meyer distribution given by equation (54). 
This assumption increases in validity with the Mach number. At a 
Mach number of 1, the flow pattern derived by Guderley (refer- 
ence 40) indicates the approximate validity of such an assumption, 
except in a small region near the beginning of the supersonic 



40 


HA.CA TN 2130 


velocity distribution, that is relatively more substantially influ- 
enced by compression waves from the sonic boundary. The preceding 
assumption is therefore probably a much better one in the present 
speed range than in that of the preceding section. 


As in the preceding section, by assuming a Prandtl-Meyer dis- 
tribution the problem of calculating the supersonic part of velocity 
distributions in the upper transonic range is again reduced to deter- 
mining one point of this distribution at each free- stream Mach num- 
ber Mq in the range, which is necessary to evaluate the constant 

Q & in equation (54). This calculation is performed along the same 

lines as in the preceding section. The continuity condition is first 
set up for the flow across a section AB in the y-direction (fig. 9), 
where A is a point on the shock, to be specified presently, at a 
distance y w from the x-axis and B is the laterally opposite point 

on the airfoil at the same chordwise location x y . Thus 

Y = (pv-l)dy (72) 

J Y 


The irrotationality condition (equation (5)) applies also in the 
present speed range, inasmuch as fluid rotation introduced by var- 
iable entropy increase across a curved shock in the small-perturbation 
transonic case, being of the order of the third power of the Mach num- 
ber increment across the shock, is negligible. Substituting equa- 
tion (5) into equation (72) and making the small-perturbation, sub- 
stitution of equation (10) for the density yields 


Y 



iv 


(73) 


where v v is the’ velocity increment on the downstream side of the 
shock at A and A is the velocity increment at B. 

A streamline- curvature function is now assumed as 


C 


= c a 



(74) 
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where s is a parameter to be determined as a function of x w and 

K. Equation (74) satisfies the boundary condition at the airfoil 
and in addition states that the curvature is zero at point A. It is 
desirable that the streamline curvature should be zero at point A 
because if it were not, then the determination of the curvature at 
point A would greatly complicate the problem. Now a point A exists 
on a curved shock at which the streamline curvature just downstream 
of the shock is zero. It is called Crocco’s point (references 38, 40, 
41, 45, and 46) and is defined by the condition 
(dv/c)0) =08 (appendix C). Instead of choosing 

' ' 'along streamline 

point A as Crocco's point, however, it will be chosen as the sonic 
point, that is, as that point at which the local Mach number just 
'downstream of the shock is unity. The reason for this choice is 
the following: Because of the additional degree of freedom pre- 

sented by the position of the shock relative to the airfoil, that is, 
determination of the distance y w , one more condition in addition 

to those already given is necessary' to make the problem determinate. 

If point A is the sonic point, a very simple such condition is 
obtained, because at the. sonic point, (d0/dv) & ^ oll g ^ = 0 (appen- 
dix C). Hence, to the first order in the change of velocity between A 
and B (that is, in the small-perturbation case) the stream direc- 
tion 9 V at A will equal the airfoil inclination 0 Q at B. The 

violation of equation (74), because the curvature will not be zero 
at A, is probably not serious inasmuch as it is known that the sonic 
point and the point of maximum streamline deflection are quite close 
to each other on the shock polar and Crocco’s point is between these 
two points (appendix C). Hence the streamline curvature should be 
close to zero at the sonic point A in the physical plane. It is 
noted that choice of a point with a particular property on the head 
shock to determine the path of integration in equation (72) makes 
the chordwise locatiort x w a function of the Mach number parameter K 

in the analysis. Hence the parameter s in equation (74) becomes a 
function only of K. 

Substitution of equation (74) into equation (73) and integration 
yields 


(YC a ) 


0 a = 9 w = 


^ r ^i(A-V v)+ -^v) 


(1-s) 


(2-s) F 


(3-s) 


(A-V 


( 75 ) 
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in which the aforementioned condition 0 a = 6 V is indicated in the 
subscript on the left side. 

The distance y v of the shock point A from the airfoil point B 

is obtained by integration of the irrotationality condition (equa- 
tion (5)) with the curvature function (equation (74)) as 

" T " y w “ c a (l-s) 


In order to satisfy the condition of zero curvature at V = by 

equation (74) and for convergence of the integrals leading to equa- 
tions (75) and (76) , it is necessary to restrict the range of values 
for the parameter s to 


0 < s < 1 (77) 

The distance d of the normal element E of the shock from the 
leading edge 0 of the airfoil can be expressed as follows: The 
angle of inclination a of the tangent AC to the shock at the sonic 
point A is given by (appendix C) 


cos a = VtV2 (78) 

This inclination is somewhat steeper than that for a Mach wave, for 
which cos a = as is to be expected. The distance CD is thus 

given by 


CD = y v cot <7 J5 y w COS a (79) 

Introducing a factor A to allow for the reduction CE due to cur- 
vature of the shock between A and E therefore gives, using 
equation (76), 

r_ AMp./2 (A.-l> w ) 

d +2v = Afn/2y v = ; (80) 

C a (l-s) 

The preceding equations will now be applied to the thin biconvex air- 
foil. For this contour, conditions at the chordwise location 

a na lyzed here are given by equation (59) as 
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(81) 


(82) 


It is convenient to proceed with the equations expressed explicitly 
in transonic- similarity form. At the sonic point behind a shock 
(appendix C), 


^ = 

- I 

(83a) 


2 

re 


775 

(83b) 


Substituting equation (81) into (80), equation (82) into (75), noting 
that the second term on the right side in equation (75) vanishes 
because' of equation (83a), and expressing equations (75), (81), and 
(80) in simi l arity form with the aid of equations (83a) and (83b) 
and the definitions 


|i/2(ir) 2 / 3 = -K/2 

(84) 

l , r^ 5 a/t 2 / 3 

(85) 


yields, respectively. 



The par ame ter s in equation (86) is next determined as a function 
of To begin with, as — >0, s must— >1, and in such a way 
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that K]_ /(l-s ) — > constant, for if not, then at E]_ = 0, Ly. < 0 

(by equation (86)), which is incorrect in view of the direction of . 
curvature of the flow. It will he assumed that the variation of s 
with is given, not only near E^ = 0, hut throughout the entire 

upper transonic range 0 < E^ < 3/2 hy 

E 2 

— — = B (a constant) (89) 

-L - ■ S 


The value of this constant can he determined as follows: Equa- 

tion (86) is similar in structure to the corresponding subsonic 
equation (56). It was assumed in the preceding section that equa- 
tion (56) should comply with the potential- limit condition 



at x = x^ = 0.5 throughout the lower transonic range. The sub- 
script x in equation (90) indicates that x is held constant in 
the differentiation. Now, if a unique continuous solution exists at 
Mo = 1, ' tllen equation (86) should become equivalent to equation (56) 

at Mq = 1 and hence should obey the potential-limit condition 
analogous to equation (90) at Mq fc 1 or E]_ = 0 and x = 0.5. An 

analogous condition to equation (90), rather than the same condition, 
is necessary because here x = x„ is a function of E-l and so can- 
not be held constant in differentiating L. However, inasmuch as 



and the second term on the right side is finite and, in fact, zero at 
Ei = 0 by equation (87), it is seen that the analogous potential- 

limit condition to be applied to equation (86) is 

) = 0 ' (92) 

E x = 0 
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The chordwise location at which equation (92)’ is applied is, at 
% = 0, "by equation (87), 2^ = 0.5 , which joins on continuously 

with the branch-point location for the biconvex airfoil in the lower 
transonic speed range. On an airfoil for which the 6 0 = 0 chord- 

a 

wise location differs from the branch-point chordwise location X2 

previously derived, interpolation between these two locations as a 
function of K would be necessary in order to obtain the branch- 
point locations in the lower transonic asymmetric-flow range. Apply- 
ing the condition of equation (92) to equation (86) results in, by 
equation (89), 



Inserting this value into equation (86) yields 


**w ,2 

and, by equation (93), 


(™a)e .0 

a 

L. > 


-11/3 


1/3 

= 2 ' = 1.260 


(93) 


(94) 


B -pTs ' 2 - 380 < 9S > 

The remaining constant A, occurring in equation (88), is 
evaluated by the boundary condition that the shock distance d— > 0 
as * K-^ a , where K-j^ a corresponds to the Mach number for 

attachment of the shock and is given by equations (71), (59), and 
(84) as K-j^ a = 3/2. With these values for d and and by 
equations (89) and (95), equation (88) yields 

(^+3/2) A = 0.126 (96) 

On the other hand, equations (89), (95), and (86) yield a cubic for 
Ly+3/2 at K]_=3/2, with the solution 

IV + f = 0.1318 (97) 


Hence 


A = 0.955 


(98) 
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The constant A is of" the expected order of magnitude, namely, some- 
what less than unity, although perhaps too close to uni ty to indicate 
a high degree of quantitative accuracy. 

The other boundary condition on shock distance, d — > oo as 
— >0, is already satisfied hy equation (88), which gives a 

limiting variation, using equation (89) and the preceding nume rical 
values. 


lim d — > 

K-j— >0 


ABL Z 0.716 


(99) 


The variation of s with K^, given hy equation (89), and the 

various numerical constants just derived enable such quantities as 
the longitudinal shock distance d, the lateral shock distance 
y v , the chordwise location x^., and the velocity L^. at chordwise 

location Xy. on the airfoil to be derived as functions of K-^ from 
equations (76) and (86) to (88). From L^K^), x^Kq.), which con- 

stitute the determination of a known point of the velocity distribution 
at each K-^, the supersonic part of the velocity distribution on the 
symmetric biconvex airfoil is obtained by the Prandtl-Meyer 
formula (54) . 


The various shock distances calculated as indicated are shown in 
figure 10. Some velocity distributions are shown in figure 7. The 
numerical values in figures 7 and 10 are plotted for a biconvex air- 
foil of =0.1. These values, however, were actually derived from 
the similarity parameters calculated by the equations of this section, 
which yield only li m iting values for t = 0, M Q = 1. As an indi- 
cation of the error involved in applying the results of figures 7 and 
10 to biconvex airfoils of finite thickness ratio, figure 10 indicates 
the Mach number for shock attachment to a t = 0.1 biconvex airfoil 
as 1.31; whereas the exact oblique-shock relations show it to be 1.48. 
Practical application of the results of this section therefore requires 
further investigation in the direction of the counterpart of the 
present calculation for t and p finitely different from zero. 
(Calculations by integral methods of detached- shock configurations 
relative to thick bodies have recently been made by W. E. Moeckel, 
reference 47 . ) 

Some comparisons with experiment. - Figure 7 has been cross- 
plotted in figure 11 to show the velocity increments A at the 
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2 = 0.5 and z = 0.8 chordwise locations through the entire tran- 
sonic speed range. From these velocity increments, local Mach num- 
bers M have been calculated by equation (27), assumi n g Ty = T= 1.2 

for over-all consistency vith the determination of A. These Mach 
numbers are shown as the solid lines in figure 12. The relative con- 
stancy of M at z = 0.5 in' the transonic range may be especially 
noted, (in the low subsonic and high supersonic ranges, M should 
approximately equal M Q . ) This phenomenon is of particular interest 
in view of the rather independent derivations of the maximum local 
Mach number at the potential-limit Mach number Mq^ and at sonic 

speed Mq = 1. The potential-limit derivation depended on the 

incompressible velocity distribution for the parameter r; whereas 
the sonic-speed derivation depended on a property at the sonic point 
of the shock polar (Be/dt>)y = 0. The relatively constant value of 

M at x = 0.5 throughout the transonic speed range is plausible 
from the expected tendency of the flow to maintain a maximum average 
value of pv ~ 1 along the line x = x^ , when both subsonic and 

supersonic velocities occur along this line. 

The discontinuities in the theoretical curves for x = 0.8 at 
Mq = 0.83 and at Mq = 0.94 correspond to the assumed sudden 

appearance of the terminal shock at the potential-limit Mach number 
(and the condition Pg = 0) and to the shock passing through the 

x = 0.8 location, respectively, (fig. 8). 

Pressure distributions in the lower transonic range on 6- and 
12-percent thick symmetric biconvex airfoils of 3-inch chord spanning 
a two-dimensional tunnel 20 inches high and 2 inches wide were meas- 
ured in references 3 and 4. The Beynolds number based on airfoil 
chord was of the order of one million. The results were not cor- 
rected for tunnel-wall effect. Terminal- shock locations plotted in 
figure 8 were estimated herein from the given pressure distributions 
and schlieren photographs. The local Mach numbers given in refer- 
ence 4 were presumably obtained by applying the finite-perturbation 
formulas (equations (6) and (26)) to the pressure differences actually 
measured. Accordingly, these formulas were used in obtaining velocity 
increments and local Mach numbers from the experimental results given. 
All velocity increments A, local Mach numbers M, and free-stream 
Mach numbers Mq were reduced from T = 0.06, 0.12 to T = 0.10 by 

the transonic similarity rules of equations (24), (25), and (31). 

The experimental results as thus reduced are shown in figures 8, 11, 
and 12. 
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Several points imcorrected by tbe transonic similarity rules 
are included in figures 11 and 12. It is seen that these rules do 
substantially correlate tbe different thickness-ratio results. The 
residual discrepancy in the results for x = 0.5 after correlation 
may be due to the finite difference between Mq and 1, to the 
finite-perturbation velocities involved, to differential tunnel-wall 
effects, or to differential Reynolds number effects. Also included 
in figures 11 and 12 are several points similarly reduced from the 
tests described in reference 48. These tests were made on a 6-percent 
thick, 4- inch chord symmetric biconvex airfoil in a two-dimensional 
tunnel 18 inches high and 4 inches wide. Beynolds number was of the 
order of one million. The discrepancy between the highest speed 
points of reference 48 and the lowest speed points of references 3 
and 4 is possibly due to larger tunnel-wall effects and the somewhat 
higher Beynolds number in the tests of reference 48. 

Agreement of the present theory with the experimental results 
is considered to be good at least qualitatively. The theoretical 
calculations are, as previously noted, limiting calculations for 
T = 0 and Mq = 1. As an indication of the necessity of allowing 
for finite departures from at least the transonic limit Mq = 1 in 
order to attain closer agreement with experiment, several alternative 
theoretical local Mach number curves are included in figure 12. In 
one, the function of equation (28a) was used in equation (27) 

to obtain local Mach numbers from the theoretical velocity increments 
of figure 11. In the other, the exact equation (26) was used. It 
Is seen that the differences between these alternative curves is of 
the same order as the difference between theory and experiment. The 
experimental local Mach numbers appear, nevertheless, to be definitely 
lower and higher than the theoretical values at x = 0.5 and x = 0.8, 
respectively. It remains to be seen whether experiments at higher 
Beynolds number would not better the agreement with the (essentially 
Infinite Beynolds number) theory. The changes with increase of 
Beynolds number may be expected to be in the direction for better 
agreement. 


FBESSUKE BRAG; USE OF MOMENTUM INTEGRAL 

The flow patterns obtained in the preceding sections, together 
with the momentum equation in integral form, yield an estimate of 
the pressure drag of the symmetric biconvex airfoil in the transonic 
speed range. The pressure drag of an airfoil is defined as that 
portion of the total drag that can be obtained by an integration 
around the airfoil of the component in the free- stream (x-) direc- 
tion of the normal force or pressure on the airfoil. The remaining 
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portion of the total airfoil drag, the skin friction, due to the 
tangential forces on the airfoil surface produced by viscous action 
in the boundary layer, will not be considered here. The skin fric- 
tion, for example, on a 10 -percent thick airfoil, is known to be 
small relative to the pressure drag in the transonic range. 

In the range of continuous potential flow, that is, from Mach 
number zero to the potential-limit Mach number M Q the pres- 
sure drag is zero, corresponding .to the fact that the completely 
subsonic potential-flow field bounding the airfoil provides no 
mechanism for transporting x-direction momentum away from the air- 
foil. In the range of transonic flow with shock, however., the shock 
structure provides such a mechanism, both in the direct shock wake 
of decreased total head and in the additional wake due to boundary- 
layer thickening in traversing the term in al shock. In this range, 
a nonvanishing pressure drag is therefore obtained that could be 
directly calculated by integration over the airfoil of the x-component 
of the pressure. This direct integration can be carried out here over 
the rearward part of the symmetric biconvex airfoil for which velocity 
distributions have been previously derived. The contribution of the 
• remaining forward part of the airfoil to the pressure drag can be 
obtained from knowledge of the conditions in the flow field by use of* 
the momentum integral. 


Integration of the x- component of the pressure between two chord- 
wise locations x a and X 5 yields for the pressure-drag coefficient 

c d,r' of 'tki 8 part of airfoil (upper and lower surface) 


r x b 


c d,r = 


D 


= 2 


2 p o v o c 


P tan 0 a dx 


= -4 


J % 




'*b 


Ae a dx (100) 


in the small-perturbation case. For the symmetric biconvex airfoil 
equation (100) becomes, by equation (59a), 



(101) 


Substituting the transonic small-perturbation limiting form (54) of 
the Prandtl-Meyer velocity distribution into equation (101) yields 



Ol 

o 


c 


d,r 


5/3 

473 




1/3 




In the lower transonic range, equation (102) was used to calculate the pressure drag of 
the "biconvex airfoil T = 0.1 from the midpoint £ a = 0 to the terminal- shock location 

=■ The pressure drag from £ = £^ to the trailing edge f 3 1 is zero, by the 

assumption made in the section LOWER TRANSONIC ASYMMETRIC FLOW WITH TERMINAL SHOCK that 
free-stream static pressure exists on the airfoil downstream of the terminal shock. In 
the upper transonic range, equation (102) was used to calculate the pressure drag of 
the biconvex airfoil from the chordvise location equal to that of the sonic point of the 
detached shock =■ to the trailing edge ■ 1, 


The pressure- drag contribution of the remaining forward parts of the biconvex airfoil 
in these ranges was determined by the momentum Integral. Thus, the net decrease of flux 
of x-momentum out of the region ABIE, (fig. 2) (momentum flux is here defined as statio 
pressure plus flux of kinetic momentum) gives the pressure drag on the portion of the 
airfoil from the leading edge 0 to point A. Ab shown in detail in appendix B, the 
resulting expressions for the pressure drag on the biconvex profile are as follows: 


(a) Lover transonic range; leading edge, £ 3 -1, to midpoint £ - 0 

5/3 p- 

c d,f = ^73 l3 [“ 2(3?FT + i T4^jJ (103) 

(b) Upper transonic range; leading edge, £ ® -1, to £ » £ , L » b 


T 5 / 3 It V , 

2 

(L^ + Ki) 3 

+%) 2 s 

c a,f = * L w + K i) 

3 

(4-s) 

(3-s) + 6 Il^Tj 


(104) 


KACA TN 2130 



NA.CA TN 2130 


51 


The resulting pressure drag on the symmetric biconvex airfoil 
is 


jjj c d = c d,f+ c d,r ( 105 ) 

cn 

and is shown in figure 13 as a function of free- stream Mach number. 
The finite value of c^ = 0.015 at the potential-limit Mach num- 
ber Mq = 0*.83 results from placing the terminal shock at the 

potential- limit chordwise location at this Mach number (see pre- 
ceding section) . As the terminal shock moves toward the trailing 
edge', the drag rapidly rises to a maximum value of c^ = 0.098 

at sonic speed Mq = 1. There are several points of interest about 

this value. First, it is of the same order of magnitude as that 
derived by Guderley (reference 40) by an entirely different (hodo- 
graph) method for a somewhat differently shaped 10-percent- thick 
airfoil. Second, the value of c^ f = 0.017 given at M q = 1 by 

the lower transonic formula (equation (103)) was slightly different 
from the value 0.011 obtained at Mq = 1 from the upper transonic 

formula (104). This slight difference is due to r approaching a 
different value, 2, from s( — *1) as Mq— >1. Inasmuch as this 

difference in c^ f is not great, and in view of the considerable 

differences in procedure on the two sides of sonic speed, the small 
difference in c^ f is to be regarded not as a defect but rather 

as indicating the consistency of the present approach. The occur- 
rence of a drag maximum at Mq = 1 rather than somewhere in the 

vicinity is probably due to the approximations involved in the 
present calculation, especially that which assumes that the terminal 
shock starts at the trailing edge. 

As the free- stream Mach number increases from sonic speed, the 
drag coefficient decreases because of the decreasing over-all pres- 
sure level. As the shock attaches and the field of flow becomes 
completely supersonic, the drag fin ally ap proaches the Ackeret rate 

of variation with Mach number, 1 //\|m q 2 -1. Several drag coefficients 

have been calculated in the completely supersonic range by using the 
exact Prandtl-Meyer pressure distribution with initial conditions at 
the leading edge given by the exact oblique- shock relations of refer- 
, ence 49. The difference between the first such point, c^ = 0.051 

at Mq = 1.50, at which sonic speed occurs behind the attached shock 
at the leading edge, and the point c^ = 0.065 at Mq = 1.31 by the 
present method, indicates again the approximation introduced into the 
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present results primarily by the transonic approximation jl << 1 
(that is, setting all the r values equal to 1.2") and, probably 
to a lesser extent, by the small-perturbation approximation 

A << 1. 


CONCLUDING REMARKS 

The small-perturbation form of an integral method has been used, 
with the aid of simplifying assumptions indicated by experiment or 
other theory, to calculate approximate flow patterns for a symmetric 
biconvex airfoil through the transonic speed range. The wholly sub- 
sonic symmetric flow field develops symmetric continuous locally 
supersonic regions when the subsonic free- stream Mach number exceeds 
a definite (lower critical) value. Above an equally definite (poten- 
tial limit) free-stream Mach number, no continuous symmetric solu- 
tion exists by the present method. Above the potential- limit Mach 
number, the method yields both a physically unreal symmetric solution 
containing forbidden regions bounded by limiting lines and a physi- 
cally real asymmetric solution containing an asymmetric locally super- 
sonic region. Termination of this locally supersonic region on the 
downstream side by a shock required additional hypotheses concerning 
the shock. These hypotheses were quite unrelated to the limiting- 
line phenomenon. It was found possible to join the asymmetric solu- 
tion continuously to the asymmetric solution that exists for slightly 
supersonic free-stream Mach numbers. The main results of the method, 
such as the variation of pressure drag with Mach number, relative 
constancy of local Mach number ahead of the terminal shock in the 
transonic range, and so forth, agree qualitatively with experimental 
data. More extensive experimental data are needed in the form of 
pressure distributions, wake surveys, and terminal- shock surveys on 
airfoils of simple shape under isolated airfoil conditions at high 
Beynolds numbers. The present method is probably simple and flex- 
ible enough to enable suitable alterations and extensions to be made, 
whether for more rigorous solutions of the cases treated, or to 
more general situations, such as airfoils with lift, bodies of 
revolution, and so forth. 


National Advisory Committee for Aeronautics, 
Lewis Flight Propulsion Laboratory, 

Cleveland, Ohio, Decenber 19, 1949. 
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APPENDIX A 


SMALL-PERTURBATION APPROXIMATION OF 
CONTINUITY INTEGRAL 

The exact form of the equation of continuity, equation (1), 
can he written as follows (see fig. 14): 




d Yc 




pv cos & 4y 


(Al) 


where 8 is a possible lateral displacement of a streamline far 
from the airfoil. Subtracting y 0 -Y from both sides of equa- 
tion (Al), as in the section COMPRESSIBLE POTENTIAL FLOW,' and pro- 
ceeding to the limit y 0 — >oo yields 


Y-8 


>00 


Oy 


(pv cos 0-l)dy 


(A2) 


An expression for dy in terms of v and 6 is next derived 
starting from the differential equations of motion in intrinsic 
form. These equations are (reference 35) 


Bv _ 
r — » -Cv 
on 


(A3) 


- ^(m 2 -i) 


Be 


Bn 


(A4) 


where 



(A5) 


and the subscript s denotes values taken along a streamline. 
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Consider the variation of velocity v and stream direction 0 
along a line mAk-T-ng an angle a D with the free -stream direction 
(fig. 14). If z denotes distance measured along this line, then, 
using equations (A3) to (A5), 

5v _ /5v\ /50\ 5s by bn 

bz ~ \50/ s \5sy s 5z 5n 5z 



C cos (<x o -0) - Cv sin (a o -0) 


= -Cv sin (a o -0) 


cot (a o -0) 4- 1 

L t Ws 


(AS) 


50 / 50\ 5s_ 50 5n • 

5z \dsy B 5z 5n 5z 


2 

-C cos (oq- 0) - C sin (a o -0) 


-C sin (oq-0) 


cot (Oq-0) 


Jm 2 -i) /5v 
T \ Se , 


(A7) 


Equations (AS) and (A7) yield, respectively, for variations along z, 

dv 


dz sin (o6 o -0) 


T + (zr?\ cot ( a n _e ) 


<*/« 


(as) 


dz sin (a o -0) 


d0 




(A9) 


If the angle a Q is now taken as jt/2, the line element 
dz becomes dy, and equations (AS) and (A9) become 
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dy cos 0 


&Y 


v + tan 9 



s_I 


(A10) 


dy cos 9 = 



d0 


c 

"1 „ (M 2 -l) fl 

tan 0 + ■ ■■ r 

-)] 


v V 

W s _ 


(All) 


Substituting equation (A10), the definition v = 1+u, the small- 
perturbation expressions cos 0=1- 0 2 / 2, tan 0=0, and the 
expansion equation (10) of p in powers of V into equa- 
tion (A2) yields 


Y-B 






<*>08 


dU 

C 


in which 





(A12) 


(A13) 


S (H) 8 . 

• 

The next higher order terms in the bracket of equation (A12) would 
be of order u 3 , V0 2 , pi?^0Vg a . The order of magnitude of 0Ug g 
is v. if equation (A12) is consistent with the transonic similarity 
rules, then 0 2 ~ pu 2 . The two terms involving 0 in equa- 
tion (A12) therefore can at most only affect the fraction r c (M 0 ) 
but not its limiting value 1^(1) = T = (7+l)/2. Thus the small- 
perturbation equation (4) of the continuity integral is correct to 
the order v in the subsonic range p ~ 1 and to the order v 2 
with r c in the transonic range p « 1. The small-perturbation 

form of the irrotationality condition (equation (5)) is seen from 
equation (A10) to be correct to the first order In v and 0. 

In the important special case where the airfoil has fore-and- 
aft symmetry, equations (4) and (5) hold exactly at the midchord 
location, for here, by symmetry, 0 = 0 in equations (Al) and (A10) . 
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In the neighborhood of such a location and, more generally, in the 
neighborhood of any 9=0 location on an airfoil, it may therefore 
be expected that the essential small -perturbation approximations 
are those involving the velocity increment u, as in equation (10) 
rather than those involving the flow direction . 9 . This expectation 
is the justification for the use of T c rather than F in the 
continuous -potential -flow calculations of this paper. 
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AP PKNDT7 B 


SMALL-PERTURBATION APPROXIMATION OF 
MOMENTUM INTEGRAL 


The momentum principle in integral form states that the Ye o tor 
rate of flow of momentum out of a fixed volume of space equals the 
vector-integrated pressure on the hounding surface of the volume. 
Applying this principle to the two-dimensional region ABDE (fig. 14) 
gives, for components in the x -direction, 


-y c 

Jo 

ED 


O7 0 +S 


(Po+Po^o 2 )^ 


P7n+S 


BB 


Po ay 


rtf 


(p+pv 2 cos 2 0)dy + 


nIY 

AB 


Pa d y 


'-’o 

EOA 


(Bl) 

in which the paths of integration are indicated under the respective 
integral signs. Rearranging the terms of equation (Bl) as with 
the continuity equation (Al), proceeding to the limit y Q — > cn, 
denoting by p, p, and v fractions of the respective free- 
stream quantities, and noting that P 0 v 0 2 /p 0 «= 7M Q 2 results In 



n Y 

ny G 


Y-8 at 

p a" 1 

-^-5- dY + 


+ (p T 2 cos 2 0-1) 


o TMo 2 

Y 

7Mo 2 J 


dy 


(B2) 


Subtraction of the continuity equation (A2) then yields for 
the pressure-drag coefficient c& f on the airfoil from the leading 
edge to the point A (including a factor 2 for both surfaces of the 
symmetric airfoil). 


rtf 


vco 


v d,f 


7Mq 2 


(p_-l)dy = -4 


Oo 


F* 1 - + pv cos 0(v oos 0-1)1 dy 

Jy L**? J 

(B3) 
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Substituting the definition v = 1+v, the expansion 
cos 0=1- 0^/2 + . . . , the expansions for p and p as power 
series in v obtained from equation (6), and the irrotationality 
expression (A10) into equation (B3) yields 



and Vq b is defined in equation (A14). Equation (B4) shows why 
the momentum, integral equation (B2)(or rather the derivative of 
equation (B2) with respect to Y) is inconvenient to use along 
with equation (4) in solving the velocity-distribution problem. 

The' flow-direction term 0 2 /2 in equation (B4) is of the same 
older of magnitude as the velocity tern — (p. / 2)v 2 and so cannot 
be neglected in the small-perturbation case, as was the same term 
in the continuity integral equation (A12). The other 0 term in 
equation (B4) is of order pU 3 and, as in the continuity integral, 
at worst merely modifies the function r m (Mq) but not its limit- 
ing value at sonic speed F m (l) = T « (l+7)/2 . Thus, in this paper, 
of the three first integrals of the differential equations of 
motion - that is, conservation of energy, mass, and momentum - the 
first two (the first as Bernoulli’s equation) were used in calcu- 
lating the basic flow pattern, and the x-component of the third 
integral, on the basis of the calculated flow pattern, was used 
to determine a portion of the pressure drag. The lift and the 
moment of the forces acting on different portions of the airfoil 
could similarly be calculated from conditions in the flow field 
by the y-component momentum integral and the angular momentum 
integral, respectively. 

For the calculation of the pressure drag on the forward 
portion of the biconvex airfoil in the lower transonic range, the 
forward portion extended to midchord . Hence the two terms involving 
0 in equation (B4) were considered as vanishing to a sufficient 
approximation. The resulting integration of equation (B4) using 
the curvature function (equation (33)) then yields equation (103). 

In the upper -transonic range, the forward portion of the airfoil 
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extended to the 9 = 9 V location and. 6 was held constant at 
this value in the integration of equation (B4). The term involving 
Vg s in equation (B4) vanishes in the limiting case Mq = 1 or 
r m =T, for which the calculations in this range were made. Inte- 
gration of equation (B4), using the curvature function and other 
relations pertaining to the upper transonic range (see section 
UPPER TRANSONIC FLOW WITH DETACHED SHOCE) leads to equation (104). 
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APPENDIX C 


TRANSONIC SMALL -PERTURBATION APPROXIMATION OF 
OBLIQUE -SHOCK RELATIONS 

An element of oblique-shock front, inclined at an angle o 
to a free-stream flow at Mach number M 0 (fig. 9) produces a flow 
deflection 0 given "by (reference 49) 


/ 7-r 

V +1 / 


(7-l)M 0 2 


+ 1 - cos 2 a 


sin 2 a tan (a-0) 


( 01 ) 


The corresponding velocity ratio v just downstream of the element 
is 


v 


cos a 
oos (.0-9) 


(C2) 


Elimination of a between equations (Cl) and (C2) yields the 
equation of the shock polar. A relation between the slope 
(dv/cteJah of the tangent to the shock polar, that is, the rate 
of change of velocity ratio with respect to stream direction along 
the downstream side of the shock and the similar rate of change 
along a streamline (c$v/c)0) a at the same downstream point of the 
shock, is obtained by dividing equation (B8) by equation (A9) and 
replacing a Q by a j thus. 


<*\ t + (§f) B cot (a -e) 

^ ■ oot (a . e) + S&H 


(C3) 


A similar relation for the rate of change of velocity with stream 
direction in the y-direction at a point on the downstream side of 
the shock element is given by setting a = rt/2 in(C3), 



v + tan 0 





o 


(C4) 
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For small- perturbations v- 1 = v « 1, equations (Cl) and 
(C2) become 


rv = - re tan ct 


- _ti_ + cos 2 ct 
M 2 


(C5) 


Also, making the transonic approximation cos a ~ ^M Q 2 -l = « 1 

and using equation (27) with T M » T yields, in place of equations 
(C5), (C4), and (C3), 


Te t= (j. oos a - 


cos 3 a 


(C6) 


Tu = -re/cos ct = - p + cos 2 a 


(C7) 


(bv\ 1 + cos CT (bv/be) B 

v^/sh cos ct + (ir+2Fi) (du/Be) s 

Equation (C7) can he differentiated to yield 


KiL ■ - 2 coe ° 

(dv\ 2 COS CT 

\d®/gh p -3 cos 2 CT 

Equations (C8), (C9), and (CIO) combine to giro 



(5 cos 2 CT*i-0 
cos o (7 cos 2 a — 3 jT) 



p 

(2 cos ct-p)( 5 cos^a-p) 
cos CT (7 cos 2 CT -3p) 


2 „ ~\ 


(C8) 


(C9) 


(CIO) 


(Cll) 


(C12) 


(C13) 
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Tlie ratio of the local curvature Cg^ = - (So/Bzjg^ of a shock 
front to the curvature of the streamline just downstream, of the 
shock element is obtained by dividing equation (A8) by (da)^, 
setting a Q = a, and going to the transonic small -perturbation 
limit. Thus, 

0 ^ _ J±£±_cos^o)_ (ou) 

t cos o (3p-7 cos 2 a) 

The conditions at various points on the shock polar in the small- 
perturbation transonic range can now be determined by the preceding 
equations. Crocco's point, for example, is defined by (Sp/S0) a =«>. 
Equations (C12) and (C14) show that this condition implies zero 
streamline curvature if the shock front itself has finite curvature, 
which is usually the case for a detached shock wave. Crocco's con- 
dition yields, by equation (C12), 

cos 2 a *= 3/7 p (C15) 

for the inclination of the shock front in the physical plane at 
Crocco’s point. Using equation (C15), the various flow quantities 
V , 0, and the local Mach number M just downstream of the shock 

element can be determined from equations (C6), (C7), and (27) 

with r M «= r. 

The derivative (Be/du) y (equation (C13)) is infinite at 
Crocco’s point, zero at the sonic point (defined by M «= 1), and 
zero also at a subsonic point of the shock polar defined by 
cos a ■=> l/5 jT. At this last point, the downstream streamline 
starts out at constant pressure. Inasmuch as here (ck>/c)0) a = 0 by 
equation (C12). 

The properties of various points on the shock polar in the 
small-p>erturbation transonic limit, determined in the manner indi- 
cated, are listed in similarity form in table H. Equations (64) 
and (68) to (70) are obtained from the values of Tu /jl given in 
the table. The velocity increment v in the table is (vgAl) - 1, 
and hence equals A- 1 -A 2 in the section LOWER TRANSONIC ASYMMETRIC 
FLOW. WITH TERMINAL SHOCK. S imi larly, fT in table I is -p+2rAp 
in the same section. 

Further discussion of the shock polar in the transonic small- 
perturbation limit is given by Guderley (reference 38), who shows, 
in particular, that from Crocco's point to the point of maximum 
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def lection on the shock polar, the oblique shock in the physical 
plane is attached to the vertex of a wedge with curved sides and 
has infinite curvature at the vertex (Cg^/C = co); reference 38 
also shows that the point of maximum deflection on the shock polar 
corresponds to the point of shock detachment from the vertex, 
which had been doubted by Crocco (reference 46) and more recently 
by Thomas (reference 50 ) . 


TABLE II 



Defini- 

tion 

SI 

cos 2 a 

cos a 

Pu 

V- 

re 

r 1 / 3 ^ 

H 

n 


iI3/2 

mm 

V- 

Mach 

wave 

•V a 0 

-l 

1 

1 

0 

0- 

0 

1 

Sonic 

point 

M a 1 

0 

1/2 

0.707 

-0.500 

0.354 

-1 

0 

Crocco’ s 
point 

(E = ° 

CO 

3/7 

0.655 

-0.571 

0.374 

m 

-0.143 

Maximum 

deflec- 

tion 

o 

II 

01 

(plso. 

-l/W 

1/3. 

0.577 

-0.667 

0.385 

-1.260 

-0.300 

Constant 

pressure 

(¥) "0 

0 

l/5 

0.447 

-0.800 

0.358 

-1.587 

-0.600 


Normal 

shock 
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APPENDIX D 


TERMINAL -SHOCK CONDITION IN LOWER 
TRANSONIC RANGE 

The condition that for finite perturbation of the flow the 
terminal shock in the lower transonic range be such that its back 
pressure is as close as it can get to free-stream pressure may be 
regarded as a rough statement, for a specific purpose, of the 
theorem that the equations of motion can be derived from the vari- 
ational principle 


5 JJ (Po-P) dr dy = 0 • (Dl) 

where p is local pressure at a point (x,y) and p Q is the free- 
stream pressure at infinity. This theorem has been established for 
nonviscous flow fields containing curved shocks (references 16, 43, 
and 61 and unpublished results of Chi -Teh Wang). It iB assumed here 
to hold at high Reynolds number in the lower transonic range More- 
over, the extremum value of the integral in equation (Dl) is assumed 
to be the same as for irrotational flow, namely, a minimum. 

Consider a stream tube of flow Just outside the boundary 
layer (fig. 15). In this stream tube, the flow accelerates isen- 
tropically from free-stream pressure and free-stream Mach number 
at station 0 to a supersonic condition at station 1, reverts 
through a normal shock to subsonic conditions at station 2, and 
returns ultimately to free-stream pressure, but not to free-stream 
velocity, at station 3. The flow from station 2 to station 3 need 
not be assumed isentropic. * If the normal shook is assumed to occur 
in a diverging portion of the stream tube, as is indicated by 
stability considerations (reference 52), then the maximum, velocity, 
or minimum pressure, on the airfoil occurs just ahead of the shock 
at station 1. If the shock does not form in a diverging portion 
of the stream tube, then the pressure coefficient to be derived is 
that just ahead of the shock, not the minimum in the stream tube. 

The small-perturbation assumption is not made. 

The shock back pressure P 2 is now expressed as a fraction 
of free-stream pressure p 0 . By the isentropic relation, the ratio 
of the pressure Just ahead of the shock p^ to the free -stream 
pressure p 0 is given by 
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(D2) 


Tiie ideal pressure ratio across the shock is given by the Bahkine- 
Hugoniot relation 



_2 7 
7+1 



7~1 

7+1 


(D3) 


As in the section on lower transonic asymmetric flow, to allow for 
the possibility that the pressure rise across the actual shock 
region of finite extent is different, for whatever reason, from 
the pressure rise across an ideal normal shock, a shock-pressure- 
drop coefficient is introduced: (The desirability of considering 

the influence of the shock-boundary-layer interaction in this con- 
nection was pointed out to the author by A. Silverstein. ) 


c 


s 


Ap s A p 8 
bi v i 2 p i M i 2 


(D4) 


such that the actual pressure ratio across the shock is given by 



(D5) 


Equations (D2) to (D5) combine to give for the actual shook back 
pressure ratio 



For given Mq and c s , this expression is found to have a maxi- 
mum at a value of M-j_ given by 
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Mn 


7+3 - c g (7+I) 


2 - o s (7+1)/2 

The maxim um shook back pressure ratio is accordingly 


(D7) 


7 

^1+~Mq 2 ) (7+1-Cg ) 


'max 


4 -Co (7+1) 


7 

7-1 


(7+i)(r+i-c s ) 


(358) 


It is now inferred, from the condition that shock back pressure be 
as close as possible to free-stream pressure, that if (pg/p 0 ) mflT <1, 

then the shock Mach number Mq is given by equation (D7), whereas 
if (P2/p 0 ) inny > 1, the shock Mach number Mq is given by equa- 
tion (D6) with p 2 /p 0 = 1. 


The ideal -shock case c a = 0 is discussed first. Here the 
free-stream Mach number Mq, at which (pg/pQ^ = 1, is 0.78. 
Hence, by equation (D7), the shock Mach number Mq for free-stream 
Mach numbers Mq < 0.78 is given by 



(7= 1.4) 


(359) 


For free-stream Mach numbers in the range 0.78 < Mq < 1, equa- 
tion (D6) with Pg/p 0 = 1 and- c 8 = 0 gives a relation between 
Mq and Mo. The preceding relations between Mq and Mq can be 
converted into relations between the shock pressure coefficient 


y p o 

jh > T ° 2 



(DIO) 


and Mo by equation (D2). The result is shown in figure 16. 

Shown for comparison are the pressure -coefficient curves for sonic 
velocity (Mq » 1 in equation (D2)) and for perfect vacuum (Mq.= « 
in equation (D2 ) ) . 
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Included in figure 16 are some experimental data given in 
references 44 and 53. The data of reference 53 consist of the 
minimum pressures measured on the wings of various airplanes in 
both low-speed supercritical flight (stalled pull-ups) and in high- 
speed supercritical flight (dive pull-outs). Of the six circle 
data points in figure 16, the three points at free -stream. Mach 
numbers of 0.42, 0.445, and 0.645 represent mlnlmums of all the 
pressures measured on the wings of the P-47C-1, the XP-51, and the 
SB2C-1 airplanes, respectively, in stalled pull-ups. The cluster 
of three points at approximately Mp = 0.73 represents minimum 
pressures selected in reference 53 from the pressure distributions 
measured on the wings of the XP-51 airplane in successive dive 
pull-outs . A further dive pull-out point for the XF2A-2 airplane 
coincides with the middle one of this cluster of three points . 

The data of reference 44 were obtained on a BACA 4412 airfoil 
spanning the Langley 24-inch high-speed closed circular tunnel. 

The variation of minimum pressure on the upper surface with Mach 
number at an angle of attack of the airfoil of 1°52.5* was taken 
as representative of their results. Further wind-tunnel data of 
the same type are given in reference 24. 

Comparison of the experimental data with the theoretical limit 
curve shows that the limit pressures indicated by the theory are, 
in fact, approached or reached but not exceeded. The leveling off 
of the pressure distribution of reference 44 somewhat below a local 
Mach number of 1.5 may be due to low Reynolds number, wind-tunnel 
effect, or both. Thus the concept of a limit curve derived from a 
shock condition independent of Eeynolds number and airfoil shape 
appears to have seme significance, at least in the limit of large 
Eeynolds numbers (single normal terminal shock), isolated airfoils, 
and for an airfoil shape that permits high local Mach number super- 
sonic regions to develop at low free -stream Mach numbers. 

Consider now the effect of shock pressure-drop coefficient c s . 

As noted in the section L0WEB TRANSONIC ASYMMETRIC FLOW WITH 
TERMINAL SHOCK, values of c s in the approximate range 0.1 to 0.5 
occur experimentally. Assume for definiteness a value c ? = 0.33. 

The effect of this value of c s in equations (D6) to (D8) is to 
Increase the shock Mach number Mq yielding maximum pg/pp to 1.5 
and to shift the free-stream Mach number M 0 at which (pg/pp)^^ *= 1 

to Mp = 1. Hence, the portion ABCD of the c s «= 0 limit curve 
(fig. 16) would be eliminated and replaced by the curve BE, slightly 
shifted to correspond to Mq = 1.5 instead of 1.483. The intro- 
duction of a shock pressure -drop coefficient c s thus yields prac- 
tically the same value of limit shock Mach number, and hence mini- 
mum pressure, as when c s = 0 in the range Mp <0.78. It may be 
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noted that the value o s = 0.33 yields the pressure ratio 
P 2 /P 0 ** 0*527 in equation (D6), that is, the pressure ratio cor- 
responding to locally sonic velocity in isentropic flow, which 
provides an interpretation' for the statement occasionally made 
that sonic velocity occurs "behind the terminal shock. 

The shock condition as thus far developed applies to finite 
and, in fact, quite large velocity perturbations and -where 1-M 0 . 
can, and possibly should, be of. the order 1. The transition to 
the transonic small-perturbation case is equivalent to replac ing 
the branch AB of the limit curve (fig. 16) (or its modification 
for c B =/ 0) by the tangent to AB at the point A. This fact is 
seen from equation (67), which is the limiting form of equa- 
tion (D6) for p. = 1 -Mq^ — > 0 and Mq^-1 — >0. For constant c a 
near Mq = 1 and by the shock back pressure condition Pg = 0 in 
this range, equation (67) yields -2A ] _(= P x ) as a linear function 
of 1 -Mq, tangent to equation (D6) at Mq «= 1. 

As remarked in the section LOWER TRANSONIC ASYMMETRIC FLOW WITH 
TERMINAL SHOCK, application of the shock condition pg /p Q = 1 or 
?2 = 0 In equation (67) implies approximately zero pressure change 
along the airfoil downstream of the shock. If this condition 
exists regardless of airfoil shape, a severe separation of the 
boundary layer must occur in traversing the shock. Similarly, for 
finite perturbations at lower free -stream Mach numbers, the present 
analysis applies only if the boundary layer can so vary In thick- 
ness as to provide sufficient flexibility in the airfoil boundary 
condition for the terminal shock to be free to stabilize at dif- 
ferent local Mach numbers along the airfoil. Only then does the 
possibility arise of choosing one of these local Mach numbers In 
accordance with equation (Dl) for shock location, independently of 
airfoil shape . 

With such freedom of shock location, equation (Dl) also indi- 
cates that the terminal shock will form in a diverging portion of 
stream tube, as this provides less variation of static pressure 
than if the shock formed (at the same local Mach n umb er) in a con- 
verging portion of the stream tube, which would necessarily be pre- 
ceded by a diverging part. Where substantial pressure change does 
occur downstream of the shock, the terminal- shock condition given 
herein requires modification; for example, in the small-perturbation 
case the . Pg term in equation (67) would have to be 'retained. 

Finally, it is of interest to note that th9 shock Mach number 
given by equation (D9) is also that for which the streamline curva- 
ture vanishes just downstream of a- normal -shock element of finite 
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curvature (references 18 and 43). Ihis property is derived by 
finite-perturbation considerations analogous to those of appendix C 
thus no close connection with the present derivation of equa- 
tion (D9) appears to exist. Inasmuch as a separating boundary 
layer traversing a shock should, however, tend to produce zero 
streamline curvature downstream of the shock element, this property 
of equation (D9) may also help to explain terminal-shock stabili- ' 
zation at approximately Mq = 1.5. , 
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APPEKDU E 
SYMBOLS 

The more important symbols used in this report are as follows: 
streamline curvature, - 

chord length, taken as unity 

D(drag)/jjr p 0 v 0 2 c, section drag coefficient 
X 2 

Ap a /— p-jv-^ , shock-pressure -drop coefficient 
2 

drag 

distance of normal element of detached shock from leading edge 

-2Ki = (1-^ 2 )/(tD 2 / 3 , transonic similarity parameter, 
(equations (12c) and (84)) 

pL/3 A/T 2/3, 

velocity parameter, (equations (58b) and (85)) 

local Mach number, ratio of local velocity to local velocity 
of sound 

(p-P 0 )/|p 0 T 0 2 ^ pressure coefficient 

pressure (fraction of free -stream value) 

curvature function parameters, (equations (33) and (74)) 

1 + A, resultant velocity at airfoil (fraction of free-stream 
value) 

1 + resultant velocity (fraction of free-stream value) 

— (l+£), chordwise location, distance in free-stream direction 

u 

ordinate of airfoil 

lateral distance perpendicular to free-stream direction 
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z = v/A, (distance in appendix C) 
n Mo 2 2-*y ? 

r = 1 + ~2 Mo , (equation (llo) ) 

T c = (7+l)/2, (equation (lid)) 

Pm h Mo 2 ^1 + ^ Mo 2 ), (equation (28a)) 

7 ‘ ratio of specific heats, 1.4 

5 streamline displacement at infinity, figure 2 

9 angle "between local, stream direction and free -stream 
direction, figure 2 

X = (equation (37a)) 

V3-r; n 

(i = - U = 1-Mo 2 

p density (fraction of free -stream value) 

0 angle "between shock element and free -stream direction 

T thickness ratio, ratio of maximum thickness to chord length 

X = — ~ — ° a , (equation (37"b)) 

(3-r) 2 n 3 

~ of the order of magnitude of 
Subscripts : 
a at airfoil 

e calculated sonic point of Prandtl -Meyer distribution 

1 incompressible or low speed 

7 potential limit 
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o free -stream condition 

w sonic point of oblique shock 

12 upstream and downstream of terminal shock, respectively 
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Figure 2. - Symmetric compressible potential flow. 
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Figure 3. - Plot of equation (38) 
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(a) Kaplan section* 


Figure 5. - Effect of Mach number on maximum velocity increment of 
various T = 0.1 sections, x a 0.5. 
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(c) Elliptic section. 

Figure 5. - Concluded. Effect of Mach number on maximum velocity 
increment of various T = 0.1 sections, x s 0.5. 
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Figure 6. — Jacobian of transformation from physical to hodograph 
plane. Symmetric biconvex airfoil, T = 0.1. 
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Figure 7. - Velocity distributions on T =0.1 symmetric biconvex 
airfoil in subsonic and transonic range. 
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Figure 8. - Terminal- shock location onT=0.1 symmetric 
biconvex airfoil in lower transonic range. Limiting 
calculation for T— >0, Mo — >1. Plotted by transonic 
similarity rules for T = 0.1. 
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Figure 9. - Upper transonic or detached shock range 
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Free-stream Mach number, M 0 


Figure 10. - Detached -shock location relative to 'T = 0.1 symmetric 
biconvex airfoil in upper transonic range. Limiting calculation 
for T— »0, Hq - > 1* Plotted by transonic similarity rules for 
T = 0.1. 
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Figure 12. - Local Mach number on symmetric biconvex airfoil in ■ 
transonic range. Limiting calculation for T — >0, Mq — »»1* 

■ Plotted by transonic similarity rules for T * 0.1. 
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Figure 13. - Pressure drag on symmetric bioonvex airfoil in transonic range* Limiting 

k H II v 1 DT «4*4* fl/1 I n i' 4»w > yi «f» /w» 4 a w 4 S nvil I . i w . 1 'f •» A 1 

w| juq ax« * jiuuwuu wj wa cuia wuxw ojaui.j.cu -j,uj x uxdj x vx' » — vt x« 


lO 


NACA TN 2130 





94 


NACA TN 2130 



1255 


Terminal 

shook 




Figure 15. - Representation of stream tube in lower transonic range 
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